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ABSTRACT 

Context. We present the CFHTLS-Archive-Research Survey (CARS). It is a virtual multi-colour survey which is based on public archive images 
from the Deep and Wide components of the CFHT-Legacy-Survey (CFHTLS). Our main scientific interests in the CFHTLS Wide-part of CARS 
are optical searches for galaxy clusters from low to high redshift and their subsequent study with photometric and weak-gravitational lensing 
techniques. 

Aims. As a first step of the CARS project we present multi-colour catalogues from 37 sq. degrees of the CFHTLS-Wide component. Our aims 
are first to create astrometrically and photometrically well calibrated co-added images from publicly available CFHTLS data. Second goal 
are five-band (u'g'r'i'z') multi-band catalogues with an emphasis on reliable estimates for object colours. These are subsequently used for 
photometric redshift estimates. 

Methods. We consider all CFHTLS-Wide survey pointings which were publicly available on January 2008 and which have five-band coverage 
in u'g'r'i'z'. The data are calibrated and processed with our GaBoDS/THELI image processing pipeline. The quality of the resulting images is 
thoroughly checked against the Sloan-Digital-Sky Survey (SDSS) and already public high-end CFHTLS data products. From the co-added 
images we extract source catalogues and determine photometric redshifts using the public code Bayesian Photometric Redshifts 
(BPZ) . Fifteen of our survey fields have direct overlap with public spectra from the VIMOS VLT deep (VVDS), DEEP2 and SDSS redshift 
surveys which we use for calibration and verification of our redshift estimates. Furthermore we apply a novel technique, based on studies of 
the angular galaxy cross-correlation function, to quantify the reliability of photo-z's. 

Results. With this paper we present 37 sq. degrees of homogeneous and high quality five-colour photometric data from the CFHTLS-Wide 
survey. The median seeing of our data is better than 079 in all bands and our catalogues reach a 5cr limiting magnitude of about i' AB ss 24.5. 
Comparisons with the SDSS indicate that most of our survey fields are photometrically calibrated to an accuracy of 0.04mag or better. This 
allows us to derive photometric redshifts of homogeneous quality over the whole survey area. The accuracy of our high-confidence photo-z 
sample (10-15 galaxies per sq. arcmin) is estimated with external spectroscopic data to cr^/(i+ : ) ~ 0.04 - 0.05 up to < 24 with typically 
only 1-3% outliers. In the spirit of the Legacy Survey we make our catalogues available to the astronomical community. Our products consist 
of multi-colour catalogues and supplementary information such as image masks and JPEG files to visually inspect our catalogues. Interested 
users can obtain the data by request to the authors. 
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of France, and the University of Hawaii. This work is based in part on 
data products produced at TERAPIX and the Canadian Astronomy 
Data Centre (CADC) as part of the Canada-France-Hawaii Telescope 
Legacy Survey, a collaborative project of NRC and CNRS. 
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1. Introduction 

Being the signposts of the largest density peaks of the cos- 
mic matter distribution, clusters of galaxies are of particu- 
lar interest for cosmology. The statistical distribution of clus- 
ters as a function of mass and redshift forms one of the key 
cosmological probes. Since their dynamical or evolutionary 
timescale is not much shorter than the Hubble time, they con- 
tain a 'm emory' of the initial con ditions for structure forma- 
tion (e.g. Bor gani & Guzzoll200ll) . The population of clusters 
evolves with redshift, a nd this evolution depends on the cos- 
mological model (e.g. Ek e et al. 1996). Therefore, the red- 



shift dependence of the cluster abundance has been used as 
a cosmological test (e.g. Bahcall & Fan! 1998 ; Borgani et al.l 
19991: Schuecker et al. 2003aE). A prerequisite for these stud- 



ies are large and homogeneous cluster samples with well- 
understood selection functions. Consequently, a large variety 
of systematic searches has been performed in various parts 
of the electromagnetic spectrum. The most extensive cluster 
searches and cosmological studies were performed in X-Rays 
(see e.g. iBohringer et alJ l2000t iReiprich & Bohringei 2002; 



e.g. Postman et al. 


1 19961 lOlsenetalJI 


1999; Gladders&Yee 


2000 


; Goto et al.l 


2002t iBahcall et al. 


20031 iGladders et al. 


2007 


; Koester et al 


. 2007); see also Ga 


1 (2008J) for a concise 



review of various cluster detection algorithms in the optical. 
Each of the cluster searches relies on certain cluster properties 
such as X-Ray emission of the hot intra-cluster gas or an op- 
tical overdensity of red galaxies and may introduce systematic 
biases in the candidate list creation. Hence, a careful compar- 
ison and selection with different methods on the same area of 
the sky is essential to obtain a comprehensive understanding of 
galaxy clusters and their mass properties. 

The Wide part of the Canada-France-Hawaii-Telescope 
Legacy Survey (CFHTLS -Wide) is an optical Wide-Field- 
Imaging-Survey particularly well suited for such studies. When 
completed it will cover 170 sq. deg. in the five optical Sloan 
filters u*g'r'i'z' to a limiting magnitude of i' AB =s 24.5. The 
unique combination of area, depth and wavelength coverage 
allows the application of a variety of currently available opti- 
cal search algorith ms. For instance, th e Postman matched fil- 
ter technique (see iPostman et al. 1996) applies an overdensity 
and luminosity function filter to photometric data of a single 
band survey. It can provide high-confi dence samples in the low - 
and medium redshift range (see e.g. Olse n et al. I ll999ll200ll) . 
The Red-Cluster-Sequence algorithm scans a two-filter survey 
for the Red Sequence of elliptical galaxies and is mainly used 
for t he medium to high re dshift regime with the r and z fil- 
ters ( IGladders & Yedl2000l) . The existence of five bands in the 
CFHTLS-Wide allows us to estimate photometric redshifts and 
the application of techniques using distance information (e.g. 



Mill er et al.ll2005l) . Furthermore, one of the main goals of the 



CFHTLS-Wide are weak gravitational l ensing studies of the 
large-scale stru cture distribution (see e.g. lHoekstra et aL l2006: 



Fu et al.1 120081 for recent results). This will allow us to com- 



plement and to directly compare optical cluster searches with 
candidates from weak len sing mass recons t ructions and shear 
peak detections (see e.g. ISchneiderl Il996t lErben et al. 2000; 



Bartelmann & Schneider] 12001 ; Wittman et a 


. J200JJ2003); 


Dahle et all 120031: iHetterscheidt et al. 12005 


; Wittman et al. 


2006: SchirmeretalJl2007UDietrichetal.l 20071). To perform 



these galaxy cluster studies, we perform an extensive Archive- 
Research programme on publicly available data from the 
CFHTLS-Wide. We baptise our survey the CFHTLS-Archive- 
Research Survey (CARS in the following). 

This paper marks the first step of our science programme on 
a significant area of CARS. We describe our data handling and 
the creation of multi-colour catalogues, including a first set of 
photometric redshifts, on 37 sq. degrees of five-colour CARS 
data. 

The article is organised as follows: Sect. [2] gives a short 
overview on our current data set; a detailed description of our 
complete image data handling is given in App.lAl Sect. [3] and 
Sect. |4] summarise the multi-colour catalogue creation and the 
photometric redshift estimation together with a thorough quan- 
tification of their quality. We continue to describe our data 
products (Sect. |5]l and finish with our conclusions in Sect. [6] 

2. The data 

The current set of CARS data consists of a subset of the 
synoptic CFHTLS-Wide observations which is one of three 
independent parts of the Canada-French-Hawaii-Telescope 
Legacy Survey (CFHTLS). It is a very large, 5-year project 
designed and executed jointly by the Canadian and French 
communities. The survey started in spring 2003 and is planned 
to finish during 2008. All observations are carried out with the 
MegaPrime instrument mounted at the Canada- France-Hawaii 
Teles cope (CFHT). MegaPrime (see e. g. iBoulade et al. 
2003) is an optical multi-chip instrument with a 9 x 4 CCD 
array (2048 x 4096 pixel in each CCD; a.' 186 pixel scale; 
w 1° x 1° total field-of-view). When completed, the CFHTLS- 
Wide will cover 170 sq. deg. in four high-galactic-latitude 
patches W1-W4 of 25 to 72 square degrees through the five 
optical filters u*g'r'i'z' down to a magnitude of i' AB » 24.5. 
See http : //www . cfht . hawaii . edu/Science/CFHLS/ 
and http : //terapix . iap . fr/cplt/oldSite/Descart/ 
summarycfhtlswide .html for further information on survey 
goals and survey implementation. 

Since June 2006 CFHTLS observations are publicly re- 
leased to the astronomical community via the Canadian 
Astronomy Data Centre (CADCfl At the time of writing raw 
and Elixir preprocessed images (see below), together with 
auxiliary meta-data, can be obtained 13 months after observa- 
tions. 

For the current work we consider all Elixir processed 
CFHTLS-Wide fields with observations in all five optical 
colours u*g'r'i'z' which were publicly available on 18/01/2008, 
i. e. observed until 18/12/2006. In total, the current CARS set 
consists of 37 sq. degrees split between the three CFHTLS- 
Wide patches Wl (21 sq. degrees), W3 (five sq. degrees) and 
W4 (eleven sq. degrees). The area around the defined patch 
centres (Wl: Ra=02: 18:00, Dec=-07:00:00, W3: Ra=14: 17:54, 



1 see |http : //wwwl . cadc-ccda . hia- iha . nrc- cnrc . gc . ca/| 
| cade/ 
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Fig. 1. Layouts of the current three CARS: The CARS data of this work are split up in the three CFHTLS-Wide patches Wl (21 sq. 
degrees; patch centre: Ra=02:18:00, Dec=-07:00:00), W3 (five sq. degrees; patch centre: Ra=14:17:54, Dec=+54:30:31) and W4 
(11 sq. degrees; patch centre: Ra=22:13:18, Dec=+01:19:00). In areas covered by thick lines spectra from various surveys are 
publicly available for photo-z calibration and verification (see text for details). 



Dec=+54:30:31 and W4: Ra=22:13:18, Dec=+01: 19:00) are 
covered on a regular grid with pointed observations. Names 
of individual CARS pointings are constructed like Wlmlp2 (read 
"Wl minus 1 plus 2"; see also Fig. [TJ. They indicate the patch 
and the separation (approximately in degrees) from the patch 
centre, e. g. pointing Wlmlp2 is about one degree west and 
two degrees north from the Wl centre. The overlap of adja- 
cent pointings is about 3'.0 in Ra and 6'.0 in Dec. The exact 
layout of the CARS survey fields is shown in Fig. Q] All three 
patches are covered by spectroscopic surveys which allow us 
to calibrate and to verify photometric redshift estimates later 
on; Wl and W4 overlap with the VIMOS VLT Deep Survey 



[VVDS; see iLe Fevre et aD d2005l) : iGarilli et al.1 d2008hl. W 3 
with the DEEP2 galaxy redshift survey (see lDavis et al. 2007 ). 
Moreover, CARS data from patches W3 and W4 have com- 
plete overlap with the Sloan-Digita l-Sky-Survey (SDSS; see 
e. g. lAdelman-McCarthy et al.ll2007t) . From patch Wl only the 
southern pointings Wlp3m0, Wlp4m0 and Wlplml have SDSS 
overlap. 

The Elixir data prepr ocessing performed at CFHT (see 



Mag nier & Cuillandrd 12004) includes removal of the instru- 



mental signature from raw data [bias/dark subtraction; flat- 
fielding; fringe correction in i' and z' data] and absolute pho- 
tometric calibration [determination of zeropoints, colour terms 
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and extinction coefficients, corrections for scattered light ef- 
fects which lead to significant inhomoge neous photometric ze 
ropoints across the CCD mosaic (s ee e.g. lManfroid et al.ll2001 



Koc h et al.l 120 04: Regnault 2007)]. The data is accompanied 
with comprehensive information on the observing conditions 
(seeing, sky-transparency, sky -background level) for each ex- 



posure 

After downloading all data from CADC and rejecting ex- 
posures with a problematic CFHT quality assessment we fur- 
ther process the data on a pointing/colour basis with our 
GaBoDS/THELI pipeline to produce deep co-added images for 
scientific exploitation. Our algorithms and s oftware modules to 
process multi-chip cameras are described in lErben et al. I (l2005l) 
and most of the details do not need to be repeated here. For 
the interested reader we give in App.lAla thorough description 
of the CARS data handling, data peculiarities and the pipeline 
upgrades/extensions necessary to smoothly and automatically 
process MegaPrime data. In addition, a comprehensive assess- 
ment of the astrometric and photometric quality of our data, 
together with a comparison to previous releases of CFHTLS 
data can be found there. We conclude that the CARS data set 
is accurately astrometrically and photometrically calibrated for 
multi-colour photometric and lensing studies. 

In the following we give a very brief summary of the most 
important CARS data characteristics: The first products of the 
THELI processing are 185 co-added science images accompa- 
nied by weight ma ps which character ise their noise properties 

for a discussion on the 



(see e.g. Sect. 6 of Erben et al.112005 



role of weight images in the object detection process). Our im- 
age stacking procedure first automatically identifies image de- 
fects (hot and cold pixel, cosmic ray hits and satellite tracks) 
in the individual frames and assigns them zero weight in the 
co-addition process. The stacking itself is a statistically opti- 
mal, weighted mean co-addition taking into account sky back- 
ground variations and photometric zeropoin t variations in the 



individual frames (see Sects. 6 and 7 from Erben et alj [2005 



and App.lAlfor further details). Identifying and masking image 
defects in the individual images before co-addition allows us to 
obtain clean stacked images also if only very few input images 
are contributing. This is essential for the processing of CARS 
r'-band data where most pointings are covered only by two in- 
dividual exposures. As an example, we show in Fig.|2]the final 
co-added r'-band image of the field Wlp3pl. We perform nu- 
merous internal and external tests to quantify the astrometric 
and photometric properties of our data. In App. IA.3I we con- 
clude that the internal astrometric accuracy of our data, i. e. 
the accuracy with which we can align individual exposures of 
a colour and pointing, is 0703 - 0704 (l/5th of a MegaPrime 
pixel) over the whole field-of-view of MegaPrime; our ab- 
solut e astrometric fram e is given by the USN0-B1 catalogue 
(see iMonet et al.ll2.Q03D . The co-added images of the different 
colours from each pointing are aligned to sub-pixel precision 
in all cases. 

We quantify the quality of the photometric calibration of 
our data in App. IA.4l App. lA.6l and App. IA.7l First, we inves- 



tigate photometric flatness over the MegaPrime field-of-view. 
We use data from the CFHTLS-Deep survey which keeps ob- 
serving four sq. degrees over the whole five-year period of the 
CFHT-Legacy Survey. This allows us to create image stacks 
from different periods and to investigate photometric consis- 
tency. Magnitude comparisons of co-additions obtained from 
the three years 2003, 2004 and 2005 indicate uniform photo- 
metric properties with a dispersion of o m i nttU * g > r >i< ~ 0.01 - 0.02 
in u*g'r'i' and about <r inf>z / « 0.03 - 0.04 in z' ■ We attribute 
higher residuals in z' to fringe residuals in this band. 

Our absolute magnitude zeropoints are tested against pho- 
tometry in the SDSS and against previous data releases of the 
CFHTLS. Our comparison with the SDSS shows that our ab- 
solute photometric calibration agrees with Sloan to <T a b s ,g'r'i' ~ 
0.01 - 0.04 mag in g'r'i' and o- abs , z > * 0.03 - 0.05 mag for 
z' . While the calibration in these four bands seems to be unbi- 
ased, we observe, at the current stage, a systematic magnitude 
offset in u* of about 0. lmag with respect to Sloan (CARS magni- 
tudes appear fainter than Sloan). For M*-band data from spring 
to fall 2006 our analysis suggests an Elixir calibration prob- 
lem leading to offsets of 0.2-0.3 mag in u* . 

Finally, we directly compared our flux measurements with 
those of the previous CFHTLS Terapix TQQQ3 release and 
Stephen Gwyn's MegaPipe project (see Gwvn 2008). The mea- 
surements to T0003 are in very good agreement with typical 
dispersions of 0.02mag; in many cases larger scatters are ob- 
served with respect to the MegaPipe data. Private communica- 
tion with S. Gwyn suggests that several MegaPipe stacks suf- 
fer from the accidental inclusion of images obtained under un- 
favourable photometric conditions; see App. lA.7l and App. IA.8l 
for further details. 

TableQJists average properties for seeing and limiting mag- 
nitude values in our survey data. The quoted values for ex- 
posure time (we list the typical exposure time per dither, the 
number of dithered observations per colour and the total ex- 
posure time in parentheses), limiting magnitudes and seeing 
correspond to a typical field and hence give a good indica- 
tion of what can be expected from the data. The seeing values 
(SExtractor parameter FWHM_IMAGE for stellar sources) 
are the median of measured seeing values from all co-added 
science images in the corresponding filters. We note that we 
measure a seeing of l'.'O or below for all co-added CARS stacks 
except for the w*-band image of Wlp3p3 for which we obtain 
l'.'l. The limiting magnitude is defined as the 5-cr detection 



limit in a 2'.'0 aperture via m\{ m = ZP - 2.5 log(5 



X0" sky)? 



where ZP is the magnitude zeropoint, N plx is the number of 
pixels in a circle with radius 270 and cr sky the sky background 
noise variation. The actual numbers for m\{ m in Table [T] were 
obtained from the field W4p2m0. It represents a CARS pointing 
with typical properties concerning exposure times and image 
seeing. A more detailed table listing these quantities for each 
individual field can be found in App. |A.9l 



see http://www.cfht.hawaii.edu/Science/ 
CFHTLS-DATA/ exposurescatalogs . html 



The described imaging data form the basis for the subse- 
quent multi-colour catalogue creation and photo-z estimation. 
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Fig. 2. Co-added CARS data: We show the final science image of the r 1 — band observations from field Wlp3pl (left panel) and 
the accompanying weight map (right panel). The weight map shows five extended satellite tracks which were automatically 
identified and masked before image co-addition (see App. IA.2l for details). In several CARS pointings/colours individual chips 
did not contain useful data and were hence excluded from the analysis. The pointing shown suffered from this problem in the 
uppermost row. 



Table 1. Characteristics of the CARS co-added science data: The 
table shows basic average properties of our final science data 
(see text for an explanation of the columns). 



Filter 


expos, time [s] 


mum [AB mag] 


seeing ["] 


w*(m.MP9301) 


5 x 600 (3000) 


25.24 


0.87 


g'(g.MP9401) 


5 x 500 (2500) 


25.30 


0.85 


r'(r.MP9601) 


2 x 500 (1000) 


24.36 


0.79 


('(i'.AfP9701) 


7 x 615 (4305) 


24.68 


0.71 


z'(z.MP980l) 


6 x 600 (3600) 


23.20 


0.66 



3. The multi-colour catalogues 

Our procedures to create multi-colour catalogues for the 
five-band CARS data are similar to the ones presented in 
Hildebrandt et al. d2006l) where we studied Lyman-break 



galaxies in the ESO Deep Public Survey (DPS). 



3.1. Preparation and PSF equalisation 

In order to estimate unbiased colours it is necessary to measure 
object fluxes in the same physical apertures in each band, i.e. 
for a given object the same physical parts of the object should 
be measured in the different bands. Since the PSF usually varies 
from band to band we apply a convolution to degrade the seeing 
of all images of one field to the PSF size of the image with 
the worst seeing. Assuming a Gaussian PSF we first measure 



the seeing and then calculate appropriate filter functions by the 
following formula: 



O" filter.k 



-4 



o 2 -a 2 

u worst u k ■ 



(1) 



with crfii ter k being the width of the Gaussian filter for convo- 
lution of the k'th image, <x wolst being the PSF size of the im- 
age with the worst seeing, and <x k being the PSF size of the 
k'th image. By doing so we neglect the non-Gaussianity of a 
typical ground-based PSF. Nevertheless, experience with the 
DPS shows that our procedure is sufficient to estimate reliable 
colours if the seeing values in the individual colours are sub- 
arcsecond and not too different. In CARS, the seeing values for 
a pointing typically do not differ by more than 072-073 (see 
TablelAa. 



3.2. Limiting magnitudes 

The images filtered in that way are then analysed for their sky- 
background properties. For the accurate estimation of photo- 
metric redshifts it is important to have a reasonable estimate 
for the limiting magnitude at a given object position. Therefore, 
we create limiting magnitude maps from the RMS fluctuations 
of the sky-background in small parts of the image. Here we 
use lcr limiting magnitudes calculated in a circular aperture of 
2x stellar FWHM diameter. This procedure ensures that vary- 
ing depths over the field are properly taken into account in the 
colour estimation. It may well be that an object would be de- 
tected in one part of the image whereas it is undetectable in a 
different part due to the dither pattern or stray-light leading to 
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inhomogeneous depth. By assigning position-dependent limit- 
ing magnitudes to each object in all bands we can later decide 
which flux measurements are significant and which are not. 

3.3. Object detection 

The object detection i s performed with SExtractor (see 
Bertin & Arno uts 1996) in dual-image mode and we consider 
all objects having at least 5 connected pixels exceeding 2<x of 
the sky -background variation. We will base our primary sci- 
ence analyses (galaxy cluster searches and weak lensing appli- 
cations) on the /'-band data. Hence, we generate our object cat- 
alogues based on this colour rather than on a combination of all 



available colours such as a% image (see e.g. lMcCracken et al 



2003). We use the unconvolved /'-band image as the detection 
image and measure fluxes for the colour estimation and the 
photometric redshifts on the convolved frames. Colour indices 
are estimated from the differences of isophotal magnitudes tak- 
ing into account local limiting magnitudes, i.e. if a magnitude 
is measured to be fainter than the local limiting magnitude, then 
this limit is used instead of the measured magnitude to estimate 
an upper/lower bound for the colour index. 

Additionally, we also measure the total /'-band magnitudes 
on the unconvolved image so that total magnitudes in the other 
bands can in principle be calculated from those and from the 
colour indices. However, it should be noted that our approach 
to run SExtractor in dual-image-mode with the unconvolved 
/'-band image for detection will never lead to accurate total 
magnitudes in the u*g'r'z' -bands. While adding/subtracting the 
appropriate colour index to/from the total /'-band magnitude 
yields accurate total magnitudes in one of the other bands for 
bright objects without a colour gradient, it can yield strongly 
biased results in other cases. Only catalogues created in single- 
image-mode on the different bands assure a reliable estimation 
of total magnitudes. Since our emphasis here is on estimating 
colours as accurately as possible, we do not pursue this issue 
further. 

3.4. Creation of image masks 

All CARS pointings suffer from bright stars and other large- and 
small-scale astronomical features that we would like to exclude 
from the following analysis. At least we want to know the loca- 
tion and shape of those areas so that catalogues can be cleaned 
from objects in problematic areas. Of course the regions which 
need to be masked heavily depend on the science project for 
which our data are used. Our main scientific drivers for the 
CARS data are the photometric identification of galaxy-clusters 
and their subsequent investigation with photometric redshift 
and weak gravitational lensing techniques. These applications 
require the accurate determination of galaxy surface brightness 
moments to at least fourth order. Hence, we want to exclude 
all image areas in which the light distribution of faint objects 
(often confined to a very small number of image pixels) is prob- 
ably altered by other sources. Amongst such defects are: 

- Extended haloes of very bright stars 

- Diffraction spikes of stars 



- Areas around very large galaxies 

- Various kinds of image reflections 

- Tracks of asteroids 

A complete manual masking process for the large amount of 
CARS data would be a prohibitively long and man-power in- 
tensive task. We developed a software package which gener- 
ates template masks for most image features which we want 
to reject. If necessary, these automatically generated masks are 
manually optimised later. Our tools are based on the following 
ideas: 

1. Object detection algorithms such as SExtractor identify 
astronomical sources by connected areas which exceed the 
sky background noise by a certain amount. The pixel dis- 
tribution of THELI produced images of empty fields has 
mode zero after sky subtraction. Large-scale artefacts like 
stellar reflection rings lead to local deviations of the back- 
ground. By running SExtractor with a fixed background 
value of zero and a very low detection threshold of 0.6<x, 
this local variation of the background leads to a significant 
increase in the detection of spurious objects. We examine 
the SExtractor catalogue for areas of significant over- 
densities and strong gradients in the object density distribu- 
tion. Corresponding image regions are flagged as problem- 
atic. The output of the procedure is an 8-bit FITS FLAG 
image with the size of the original image (masked areas 
are '1' and unmasked areas are '0' in this image) or/and 
a saoimage/ds9 po lygon region file of masked areas. See 



Dietrich et alJ (120071) for further details on the algorithm 



and its implementation. 

Astronomical Standard Star Catalogues such as USN0-B1, 
GSC-2 or SDSS-R5 list the positions and magnitudes of 
known astronomical sources up to a magnitude of about 18. 
In the CARS data, the large majority of these objects with 
m < 16 are bright or moderately bright stars whose sur- 
roundings should be excluded from object catalogues (faint 
haloes, diffraction spikes). Moreover, stellar sources have 
well defined shapes over the complete MegaPrime field- 
of-view. The extent of the central light concentration and 
the width and the height of stellar diffraction spikes can be 
modeled as function of apparent magnitude. On the basis 
of these observations we automatically create object masks 
for stellar objects: 

- We retrieve object positions and magnitudes from 
the Standard Star Catalogues GSC-1, GSC-2 .3.2 and 
USN0-A2. We found that our selection criteria in these 
catalogues (magnitude limits, catalogue flags) result in 
slightly different source lists and hence the three sam- 
ples complement each other. Our masking is performed 
independently on all three catalogues. 

- At each catalogue position we lay down template masks 
for the central light halo and the diffraction spikes. The 
templates are scaled with (red photographic) magnitude 
to conservatively encompass the stellar areas. In addi- 
tion, for very bright stars with m < 10.35 we mask ex- 
tended stellar diffraction haloes. For MegaPrime these 
haloes have an extend of about 4'.0 depending only 
weakly on magnitude. Moreover, these haloes occur 
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with a radial offset towards the MegaPrime centre. The 
halo displacement from the stellar centre as function of 
MegaPrime position can well be described by -0.022 
times the relative position of the star with respect to the 
camera centre. 

- Finally, the masks are converted to saoimage/ds9 
polygon region files which can furthe r be processed by 
the W eightWatcher programme (see lBertin & Marrnol 
120071) to construct a FLAG JMAGE file. 
Tracks of fast moving asteroids typically show up as a se- 
ries of high S/N, lined up, short dashed and highly ellip- 
tical objects in co-added CARS images. They are present 
in the data because our strictly linear co-addition process 
does not include any pixel rejection/clipping procedure. We 
try to detect and mask them in our multi-colour object cat- 
alogues. We identify an asteroid candidate if a minimum 
number of N objects are located within 0.4 pixels from a 
line connecting any two objects within overlapping boxes 
of M x M pixels 2 . We run this algorithm for the two pa- 
rameter sets N - A; M - 100 and N - 5; M - 175 and 
merge the resulting candidate lists. This combination was 
found empirically to give good results on the CARS data 
set. For real asteroids the ellipticities of contributing ob- 
jects are usually highly aligned. As in w eak lensing theory 
(see e.g. Bartelmann & Schneider 200 lb we compute the 
two-component ellipticity 

1 -r 

( e , , e 2 ) = (cos 261, sin 26) , (2) 

1 + r 

which depends on the object axis ratio r and position angle 
6 as determined by SExtractor. The expectation value of 
both ellipticity components is zero if the ellipticities of dif- 
ferent objects are not aligned. We then compute the align- 
ment estimator 



-4 



[median(ei)] + [median^)] 



(3) 



from the ellipticities of all objects belonging to a candidate. 
We only keep asteroid candidates with A > 0.20; A > 0.24 
(first and second parameter set) in order to minimise the 
false flagging of galaxies in areas with increased object 
number density, such as galaxy clusters. These parameters 
were optimised for typical CARS seeing conditions and im- 
age depths. For the W4mQmQ field the algorithm automati- 
cally masks 30 out of 32 visually identified asteroid tracks, 
with one false positive, for an object number density of 
35/arcmin 2 . 

We note that the different algorithms are complementary to 
each other. While large-scale features such as very large galax- 
ies or image borders influence the object density, small scale 
defects from medium bright stars (diffraction spikes; outer 
extended haloes) are caught by masking known catalogue 
sources. We independently run the object density analysis on 
all five colours of a CARS pointing. However, the stellar and 
asteroid track masks are calculated for the /'-band only. The 
latter ones need some manual revision which is done on the 
basis of the /'-band image only. Hence, asteroid tracks in the 
u'g'r'z' bands are not included in our object masks. Other prob- 
lems which require manual optimisation of the image masks 



are: (1) The object density distribution analysis also masks rich 
galaxy clusters. (2) Some objects labelled as stellar source in 
the Standard Star Catalogues are galaxies. (3) For images with 
an exceptional good seeing of 076 or better the high density of 
objects leads to a significant number of false positives in the 
asteroid masking. The final masks from the individual colours 
are merged and collected in one saoimage/ds9 polygon region 
file. The masking information is also transfered to our multi- 
colour catalogues as a MASK key which allows an easy filter- 
ing of problematic sources later. Fig. [3] shows examples of our 
masking procedure. 





Fig. 3. Semi-automatic image masking: Shown is the result 
of our semi-automatic image masking for areas of the field 
W4m0m0. The polygon squares result from our object density 
analysis and the stars cover sources identified in the GSC-1, 
GSC-2.3.2 and USN0-A2 Standard Star Catalogues (upper 
panel; multiple masks around stars appear for sources identi- 
fied in various catalogues). The lower panel shows results from 
our asteroid masking procedure. 
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4. Photometric redshifts 

From the multi-colour catalogues described in the preceding 
section we estimate photometric redshifts for all objects in 
two steps. In a first pass we use available spectroscopic in- 
formation from the VVDS0 to correct/recalibrate our photo- 
metric zeropoints on a patch-wide basis. Afterwards we ob- 
tain p hoto-z estimates for our objects (see iHildebrandt et al.1 
20081) . In the following we set the minimal photometric er- 



ror to 0.1 mag in order to avoid very small purely statisti- 
cal errors for high-S/N objects and to take into account our 
estimated internal and external photometric accuracies (see 
Sect. IA.8t . Throughout this work we use MegaPrime filter re- 
sponse curves which were computed by Mathias Schultheis and 
Nicolas Regnault. They are available at http : //t erapix . 
iap . fr/forum/sh owthread .ph p?tid=13^5 

The following analysis only includes secure VVDS objects 
(marked by flags 3, 4, 23 and 24; in total these are 4463 objects 
for Wl [up to a limiting magnitude of i'^ « 24] and 9617 for 
W4 [up to i' AB « 22.5]). Here and in the following we match ob- 
jects from our source lists with those from external catalogues 
if their position agr ees to better than l'.'O. F irst we run the new 
version of Hyperz ( Bo lzonella et al. 2000jf|on 13 fields with 
overlap to the VVDsH, four of which are in Wl and n ine in W4. 
We use the CWW template set dColeman et alJI 19801) sup plied 
by Hyp erz and add two starburst templates from lKinney et al.l 
(119961) . Additionally, we fix the redshift to the spectroscopic 
redshift for every object. In this way we find the best fitting 
template at the spectroscopic redshift for every object. Hyperz 
puts out the magnitudes of the best-fit templates and enables 
us to compare these to our original estimates. We average the 
differences between the observed and the best-fit template's 
magnitudes over all objects. In this way we derive corrections 
for the zeropoints in the five bands. We only use spectra of 
galaxies with < 21.5 which have a high S/N photomet- 
ric measurement in all filter bands; these were 654 sources in 
Wl and 2158 objects in W4. The mean and the scatter of the 
corrections in the four Wl fields are Am* = -0.064 ± 0.015, 
Ag' = 0.069+0.005, Ar' = 0.027+0.019, Ai' = -0.004±0.018, 
and Az' = 0.007 + 0.007. In the nine W4 fields we find Au* = 
-0.088 + 0.011, Ag' = 0.136 ± 0.029, Ar' = 0.019 ± 0.03, 
Ai' = 0.008 + 0.023, and Az' = -0.010 + 0.014. Note that the 
photo-z code is only sensitive to colours so that the absolute 
values of the corrections in the different bands should not be 
misunderstood as pure calibration errors. Prior to the calibra- 
tion step, we did not modify the W3 and W4 u* zeropoints for 



3 spectroscopic data were obtained from http : / /cencosw . oamp . 

|fr/VVDS/| 

4 Note that there are at least two more sets of MegaPrime filter 
curves available on the WWW: On the CFHT web pages ( [http : //] 
|www.cfhtT hawaii . edu/Instruments/Filters/megapri me . 
html) and on Stephen Gwyn MegaPipe pages (http : //wwwl . 
|cadc- ccda . hia- iha . nrc-cnrc . gc . ca/megapipe/docsTj 
[filters .html) 



5 publicly available at |http: //www, ast . obs-mip. fr/users/ 
roser/hyperz/ 



identified systematic calibration problems (see App. IA.4I >. As 
all W3 and W4 fields are equally affected by it we expect that 
it is taken into account properly by our correction procedure. 
We also did not apply any galactic extinction corrections to our 
catalogues. 

For the Wl fields that do not overlap with the VVDS we use 
the zeropoint corrections from the field Wlp2p3, the one with 
the highest density of spectroscopic redshifts in the Wl region. 
Since the regions W3 and W4 show different M*-band calibration 
systematics than Wl (see Sect. IA.6l ). we correct all W3 fields and 
the two W4 fields without VVDS overlap with the values from 
W4plml, again the most densly covered field in this region. 

Then we run B ayesian Photometric Redshifts (BPZ; 
see iBenit ez 2000£| on the catalogues with the corrected pho- 
tometry using the same template set as before. The Bayesian 
approach of BPZ combines spectral template x 1 minimisation 
with a redshift/magnitude prior. The prior was calibrated from 

HDF-N observations and the Canada-France Redshift Survey 

I II 1, 

(see Lill y et al.l 1 19951) . It contains the probability of a galaxy 

having redshift z and spectral type T given its apparent mag- 
nitude m. A detailed descr iption of the code and the prior can 
be found in Benite3d2000l) . We restrict the fitting of the photo- 
z's to z < 3.9 due to the limited depth of the Wide data. The 
Bayesian redshift estimates are added to our multi-colour cat- 
alogues. Note that not all objects in our catalogues have well 
determined photometric measurements in the full u* to z' wave- 
length coverage. This can have physical reasons (e.g. high- 
redshift dropout galaxies which are fainter than the magnitude 
limit in blue passpands) or it can be connected to problems in 
the data itself (e.g. pixels without information in one of the fil- 
ter bands). Our current catalogues miss information to cleanly 
distinguish between these cases but only allow us to identify 
problematic photometry by either large photometric errors or a 
flux measurement below the formal detection limit. In all cases 
with a magnitude estimate below the limiting magnitude, or a 
magnitude error larger than lmag, we configured BPZ to treat 
the object as non-detected with a flux error equal to the lcr 
limiting magnitude. This leads to unreliable results if the large 
photometric error results e.g. from image defects and not from 
intrinsic source properties. To allow an easy rejection of such 
problematic sources each object in our catalogues obtains pho- 
tometry quality flags for all filter bands. 

The internal accuracy of the BPZ phot o-z's is des cribed by 
the ODDS parameter (see e.g. iMobasher et al . 2004) assign- 
ing a probability to the Bayesian redshift estimate by integrat- 
ing the posterior probability distribution in an interval that cor- 
responds to the 95% confidence interval for a single-peaked 
Gaussian. By rejecting the most unsecure objects with a low 
OD DS value one c an obtain much cleaner subsamples; see also 
IHildebrandt et al.l d2008l) . 

If not stated otherwise we use in quality assessments of our 
BPZ photo-z's the following subsample of our catalogue data: 

1. We reject all objects falling within an object mask (see 
Sect. El. 



6 spectroscopic data were obtained from http : / /cencosw . oamp . 



7 publicly available at|http://acs .pha. jhu. edu/~txitxo/ 
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Table 2. Statistics of the comparison between photometric 
and spectroscopic VVDS redshifts: The third column gives the 
number of uniquely matched sources between our catalogues 
and high-confidence VVDS objects (see text). The complete- 
ness value (fourth column) lists from those objects the percent- 
age with a high-confidence BPZ photo-z estimate (ODDS > 
0.9). 



i'<22.5 



1 .5 



" defined as the percentage of galaxies with (z P hot _ Zspec)/(l+Zspec) > 
0.15 

* bias and scatter of (z p hot - Zspec)/0 + Zspec) after outlier rejection 

2. We select galaxies by means of the SExtractor star- 
galaxy classifier CLASS_STAR and reject all sources with 
CLASS_STAR > 0.95. 

3. We include only objects with reliable photometry in all five 
filter bands (see above). 

4. Finally we reject all sources with ODDS < 0.9. 

Our catalogues contain in total 3.9 million galaxies outside an 
object mask (rejection steps 1 and 2) and finally 1.45 million 
sources (about 13 galaxies per sq. arcmin) with reliable BPZ 
photo-z estimates (object sample after all rejections). 

We first compare our photo-z's from Wl and W4 to spectro- 
scopic redshift s from the V VDS in a similar way as presented 
in lHildebrandt et aL ( 2008 ). Note that these spectra were previ- 
ously used to calibrate the data! Table |2]summarises the results 
indicating a homogeneous dispersion <x A ,/ (1+Z ) * 0.04 - 0.05 
and an outlier rate (defined as the percentage of galaxies with 
(Zphot - z sp ec)/(l + Zspec) > 0.15) of 1-2% up to i'^ = 24. The 
°"Az/(i+z) statistics is estimated after outliers have been rejected. 
If we perform the spectro-z vs. photo-z comparisons with the 
ODDS > 0.0 sample (but all other filters as described above) 
the dispersion is nearly unchanged while the outlier rate rises 
by a factor 3 to 8. This confirms that the ODDS parameter is a 
good selection criterion to reject outliers and to obtain samples 
of homogeneous photo-z quality up to about i'^ « 24. A plot of 
the photo-z vs. VVDS spectro-z results in the regions Wl and W4 
is shown in Fig.|4] While the figure shows an overall good per- 
formance of our photo-z estimation it reveals residual system- 



1 - 



Field 




N 


compl. 


outl. rate r£ 


Az/(l+z£ 




[AB] 


[%] 




[%] 




Wlp2p2 


22.5 


212 
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Fig. 4. Photometric vs. spectroscopic redshifts in the Wl and W4 
regions: We show in the upper panels 1349 (Wl) and 3312 (W4) 
galaxies with i'^ < 22.5, reliable VVDS flags, good photome- 
try in all five filter bands and ODDS > 0.9 (points). Triangles 
represent galaxies with < ODDS < 0.9 (117 sources in Wl 
and 170 objects in W4). Lower panels show a binned distribu- 
tion of (Zphot -Zspec) from the ODDS > 0.9 samples of the upper 
panels. 

atics. A significant tilt is present in the z p h t vs. z spec comparison 
leading to a systematic overestimation of up to 0.1 - 0.2 of the 
redshift at low z spec and to a underestimation at high z spec . The 
tilt crosses the zero axis at z = 0.5 and hence it cancels negative 
and positive contributions to statistics involving Az (see Fig.|4]i. 
The Az/(1 + z) statistics for the complete Wl (N = 1466) and 
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Table 3. Quality parameters for the comparison of CARS BPZ 
photo-z against SDSS and DEEP2 spectra. The third column 
gives the redshift range probed by the spectroscopic samples; 
see Table [2] for the meaning of the other columns. 



Field 


Survey 


z-coverage 


N 


!?[%] 


Az/(1 + z) 


Wl 


SDSS 


< z < 0.6 


44 


2.4 


0.069 ± 0.050 


W3 


SDSS 


< z < 0.6 


208 


4.8 


0.068 ± 0.042 


W4 


SDSS 


< z < 0.6 


400 


12.8 


0.080 ± 0.039 


W3 


DEEP2 


0<z< 1.5 


448 


7.3 


0.010 + 0.050 



W4 (N = 3488) samples are: Az/(1 + z) = -0.006 + 0.051 (Wl) 
and Az/(1 + z) = -0.005 + / -0.050 (W4). If we split the sample 
at z = 0.5 we obtain for z < 0.5: Az/(1 + z) = 0.03 + 0.043 
(Wl: N = 595) and Az/(1 + z) = 0.025 ± 0.037 (W4: N = 1690). 
Accordingly for 0.5 < z < 1.0: Az/(1 + z) = -0.032 ± 0.035 
(Wl: TV = 809) and Az/(1 +z) = -0.034±0.038 (Wl: TV = 1728). 

We do not try to remedy these systematics in this article 
but we will investigate it in a companion paper (Hildebrandt et 
al. in prep.). The overestimation at the low-z is mainly caused 
by the redshift prior in BPZ. It seems to give too little prob- 
ability to the low-z population in the CARS data. A modifica- 
tion of the original prior in this sense removes the observed 
bias for < z < 0.5. The high-z underestimation of our 
redshifts can be corr ected by a recalib r ation of the original 
Coleman et ail d 1980h and Kinnev et a il dl996h template sets; 



see also iFeldmann et al.l (120061) . In the following we further 
check the consistency and quality of our current photo-z sam- 
ple. 

4.1. Internal and external quality checks on our 
photo-z, sample 

Besides with the VVDS the CARS patches overlap with public 
spectra fr om the SDS jg| (W1 -W4) and the DEEPif] redshift sur- 
vey (W3; Davis et al. 2007b . Hence, we can test our photo-zs 
which were partly calibrated and verified against the VVDS 
with an independent set of spectroscopic data. The compar- 
isons are shown in Fig.|5]and Fig.|6l formal quality parameters 
are listed in Table [3] We observe exactly the same system- 
atics identified in the comparisons with the VVDS spectra: A 
systematic tilt with overestimates of about 0.05 - 0.15 in the 
low-z regime and a reverse trend for z > 0.5. Our tests in- 
dicate that the photo-z quality and remaining systematics for 
the current CARS data set are comparable in the mean for all 
fields; regardless whether the galaxies profited directly from 
a previous calibration with spectra or whether we transfered 
corrections obtained with a galaxy subset to other fields. The 
trends we see with our BPZ photometric redshifts also show 
up in a comparison with previously obtainedphoto-z estimates 
with the photometri c redshift code Le Pharq^jon the CFHTLS - 
Deep field D 1 (see lllbert et al. I l2006l) . We show a direct com- 



spectroscopic data were obtained from |http : //cas ■ sdss ■ org . 
|astro/en/tools/search/SQS . asp] 

9 spectroscopic data were obtained from http : / /deep . berkley . 

|edu/DR3| 

10 see http : / /www . oamp . fr/people/ arnouts/LE.PHARE . html 
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Fig. 5. Comparison of CARS BPZ photo-z against SDSS spectra 
(0 < z < 0.6): The plot shows N = 44 objects for Wl, N = 208 
for W3 and N = 400 for W4. Note the different ordinate in the 
lower right panel ! 
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Fig. 6. Comparison of CARS BPZ photo-z against DEEP2 spec- 
tra for < z < 1.5: shown are N = 448 common objects be- 
tween the DEEP2 redshift survey and the CARS fields W3mlm2 
and W3mlm3. 



parison of the Dl and our Wlp2p3 photometric redshift esti- 
mates in Fig. |7] The figure confirms an overall very good agree- 
ment and a systematic tilt for < z < 1 in our estimates; 
our high-confidence BPZ photo- z sample w i th < 24 has 



24558 common objects with the Ilber t et alj d2006l) catalogue. 
The latter was filtered for (z UP i<7- - Zinfio-)/(l + ZDbert) < 0.25. 
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Fig. 7. Compariso n of our B PZ ph oto-z estimates on Wlp2p3 
with those from lllbert etail d2006l) on the CFHTLS -Deep Dl 
field. The comparison includes 24558 common objects with a 
pre-filtering for high-confidence sources in both catalogues and 
i' AB < 24 (see text for details). 853 objects (3.89%) lie outside 
the plotting region of -0.7 < z p h t - Ziibert < 0.7. We show one 
point out of five for clarity of the plot; contours indicate areas 
of 0.8, 0.5, 0.25 and 0.05 times the peak-value of the point- 
density distribution. 



For the complete common sample with < 24 we measure 
Az = (zphot - zm*rt) = -0-02 ±0.11. 

Finally we perform two internal consistency checks on our 
estimates. The first one is a comparison of independent esti- 
mates from overlap objects in different CARS pointings (see 
Fig. [T]i. In Fig. [8] we show on a patch basis Az p hot = (z P hoti - 
Zphot2) for all overlap sources. The means and scatters of this 
quantity for individual patches are: Az p hot = 0.0002 + 0.0772 
(Wl; 24329 objects), Az pho t = -0.00005 + 0.0767 (W3; 621 1 ob- 
jects) and Az phot = -0.0006 + 0.0723 (W4; 16849 objects). The 
plot and the numbers demonstrate a homogeneous and robust 
redshift estimation over the whole CARS area. Note that this 
comparison includes areas with and without spectroscopic cal- 
ibration overlap. The second internal consistency check are the 
redshift distributions for the three CARS patches. Fig. [9] shows 
that they are very comparable for all three CARS patches. 



Fig. 8. Consistency of our photometric redshift estimates: We 
compare independent photo-z measurements from objects ap- 
pearing in multiple CARS pointings in each patch (Wl: N = 
24329; W3: N = 6211; W4: N = 16849). Dashed lines mark 
regions with Az p hot = 0.1, 0.2. Only 1 out of 10 points is shown 
for clarity of the plots. See the text for further details. 
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4.1 .1 . Angular cross-correlation of galaxy populations 
in different photo-z redshift bins 

In the following we apply a correlation function analysis to 
further quantify the reliability and suitability of photometric 
redshift estimates. The theoretical background will be detailed 
in Benjamin et al. (in prep.). Judging photo-z quality by cal- 
ibration with external spectroscopy provides us with an over- 



Fig. 9. Normalised distributionw of our high-confidence pho- 
tometric redshift estimates for all CARS patches (Wl: N = 
205956 for 17 < < 22 and N = 487593 for 22 < i' AB < 24; 
W3: N = 52295 for 17 < i' AB < 22 and N = 119589 for 
22 < i' AB < 24; W4: N = 104417 for 17 < < 22 and 
N = 252969 for 22 < i' AB < 24); all distributions have only 
very few objects beyond redshift 2 (not shown). 
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all picture of the dispersion of our estimates and the total rate 
of catastrophic outliers. First, this quality control is limited to 
the magnitude and/or redshift range of our external compari- 
son sample and second we often need a more detailed picture 
on the photo-z (re)distribution of galaxies in redshift space. For 
instance, weak lensing tomography studies of the cosmic shear 
effect do not require a precise redshift estimate for each indi- 
vidual galaxy. However, we need to reliably separate galaxies 
into redshift bins and we require a precise understanding of our 
errors concerning redshift misidentifications an d inter-bin con- 
tamination factors (see e.g. lHuterer et al 12003) . 



To quantify the crucial error contribution due to large sys- 
tematic redshift misidentifications we investigate the angu- 
lar cross-correlation function of galaxies in different photo-z 
slices. The cross-correlation technique has already been ad- 
voca ted as a way to recon struct the source redshift distribu- 
tion (ISchneider et al.1 120061) . or, in combination with spectro- 
scopic redshift, as a way to improve the photo-z calibration (see 
Mandelbaum et alj|2007b . 



For accurate redshift estimates we expect to see a strong 
auto-correlation in individual redshift bins and, due to the ex- 
pected photo-z scattering (cr Az » 0.05 - 0.1), a lower-level 
cross-correlation signal in neighbouring slices. But one would 
not expect to see a cross-correlation signal for slices that are 
physically far apart. Contamination by catastrophic photo-z 
failures would lead to significant amplitudes of the angular 
cross-correlation function of photo-z slices that are well sepa- 
rated in redshift. In Fig.[l0]we show a correlation slice-analysis 
of the CARS photo-z estimates for objects with < 24. It 
represents a matrix plot including the angular cross-correlation 
functions of all pairs of photo-z slices. The figure illustrates the 
expected behaviour: We observe a very significant autocorre- 
lation and a decreasing cross-correlation in neighbouring bins. 
Our analysis also shows a decent signal for the highest redshift 
bins (z > 1.5) with low-z slices. This shows that high-z tails, 
which are often observed in redshift distributions derived from 
photo-zs, are populated with true low-z galaxies. We note that 
the ODDS > 0.9 filtering applied hitherto rejects large parts 
of dubious sources with a very broad or doubly peaked pho- 
tometric redshift probability distribution. To visualise the ef- 
fect of low-z high-z contamination we relaxed our filtering to 
ODDS > 0.8 in the cross-correlation analysis. 

In Fig. [10] we show separately the correlation functions of 
the faintest objects with 23 < < 24. We see that the faint 
population of our galaxies behaves exactly in the same way as 
the complete sample. This indicates that the photo-z accuracy 
does not degrade in the lower S/N regime. 

While Fig. [10] already allows us to draw important quali- 
tative conclusions we can derive formulae for the mutual con- 
tamination of redshift bins. A complete matrix description of 
the formalism, which will allow a consistent analysis of the re- 
sults in Fig. [10] will be presented elsewhere (Benjamin et al., in 
prep.). In this paper we limit the discussion to a strict pairwise 
cross-correlation analysis, i. e. we present quantitative results 
only for cases where the whole redshift sample is split in ex- 
actly two bins. 



The basic idea is to infer the degree of contamination be- 
tween two redshift slices from the measurement of the cross- 
correlation W\2 between the bins. 

The fraction of galaxies from bin 1 present in bin 2, as a 
fraction of the true number of galaxies in bin 1 (Nj), is fa. 
Likewise the fraction of galaxies in bin 1 from bin 2 is given by 
fa , which is defined such that the number of galaxies present in 
bin 1 from bin 2 is given by Njfa ■ Hence the observed number 
of galaxies in each bin (N°) can be expressed as, 

Nj(l-.fa) + Njfa, 

2KL-f2i) + Njf n . (4) 

The first term of each equation accounts for those galaxies that 
do not leave the given bin, the second term accounts for those 
interloping galaxies from the other bin. Inverting these equa- 
tions, the true number of galaxies can be expressed as a func- 
tion of the observed numbers and the fractions fa and fa, 



N° 2 



Nj(l 



N = 



N 



N° - fa(N° + N°) 

1 - /12 - fa 
N° - fa(N° + N° 2 ) 



1 - /12 - /21 

Note that N°+N 2 - Nj +Nj, thus the total number of galaxies 
is preserved, as should be the case. It is also obvious that for 
cases where /12 +/21 is unity there is a zero in the denominator. 
What is less clear, is that in these cases the numerator is also 
zero, which can be seen by plugging Eq. [4] into the numera- 
tor. In this case the system of equations is degenerate, and will 
not admit a unique solution. This should not pose a practical 
limitation since it is expected that the fractional contamination 
between bins is small, and specifically less than 0.5. 

In order to calculate how the cross-correlation function is 
changed for non-vanishing coefficients fa, fa, it is sufficient to 
consider the natural estimator of the angular correlation func - 
tion, as opposed to that presented by lLandy & Szalayl (11993b . 
The natural estimator works well at small and intermediate 
scales where edge effects are not an i ssue, provided that th ere 
is a sufficient density of points (see iKerscher et al. 1 120001 for 
a comparison of the estimators). The observed angular cross 
correlation functions are given by, 

(PiD0° g 



1 + w° = 



(R l R2)e 
(D X D 2 )° 



(6) 



(7) 



(R\Ri)e 

where {D\D\)° is the observed number of pairs separated by 
angle within bin 1, similarly (D\D2) „ is the number of pairs 
between bins 1 and 2, and (R\R2)e is the number of pairs be- 
tween objects from random fields of identical geometry. 

Considering how galaxy pairs are split between the two 
bins 1 and 2, one can show that the observed number of pairs 
depends on a combination of the true number of pairs and the 
contamination fractions: 



(DiDife = (£>iD 2 )T((l-/ I2 )(l-/ 21 )+/ 2I /i2) 
+(DiDi)J(l-/i 2 )/i2 
+(D 2 D2)J/ 2 i(l-/2i). 



(8) 
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Plugging this relation into eq. (0, and noting that the term 
(PiD 2 )°l(RxR 2 ) e must be normalised by NfNf/N°N°, where 
TYj , is the number of objects in the random samples, the fol- 
lowing equation can be derived for the observed angular cross- 
correlation function, 

(N T ) 2 

l + a>1 2 = (1 +W T)_L-/ 12 (1-/ 12 ) 

(Nj) 2 



+d +^22)77^/21(1-/21) 



+(1 +cj ,) 



N°N° 

T A _1 2_ 

N°N° 



(I-/12-/21 +2/12/21). 



(9) 



Note that the observed cross-correlation function depends on 
the unknown true number of galaxies in the bins and the un- 
known true auto-correlation function. The true galaxy number 
can be expressed in terms of the observed number of galax- 
ies and the contamination fractions via eq. ©. It is possible to 
express the true auto-correlation as functions of contamination 
fractions, the number of observed galaxies and the observed 
auto-correlation functions (Benjamin et al., in prep), 



^22 



(1-/20 2 



(i 



/l2) 2 d-/2l) 2 

f 2 
■'21 



' /12/il 



(l-/l2) 2 (l-/2l) 2 -/i 2 2/ 2 2 l 
2/2ld-/2l) 



(l-/l2)(l-/2l)+/l2/21 



(10) 



By exchanging 1 and 2 in eq. ( fTUb an equivalent expression for 
the auto-correlation of bin 2 is obtained. To finally use eqs. (0 
and dTOb we make the explicit assumption that the true cross- 
correlation between the two redshift bins is zero (o>L = 0), i.e. 
all the observed cross-correlation is due to contamination. This 
prescription allows us to use the observed correlation functions 
and number of galaxies to determine the contamination frac- 
tions, /12 and fi\, for a pair of redshift bins. 

We note that the outlined formalism cannot be trivially ex- 
tended to a multi-bin setup, since it assumes a pair of bins and 
ignores possible contamination from other redshifts. However, 
it already allows us to recover the fraction of objects that cross 
a given redshift z cut due to photometric redshift errors, and an 
analysis can be done as a function of z cut . 

We apply the pairwise analysis on our data by cutting it at 
z cut = 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 1.0; 1.5; 2.0 yielding a 
low redshift bin 0.0 < z\ < z C ut and a high redshift bin z cut < 
Z2 < 4.0. The angular auto and cross-correlation functions from 
the pair of bins are used to estimate the contamination fractions 
fi2 and /1 by fitting the observed cross-correlation with eq. (0. 
The analysis was performed with eleven equally spaced cross- 
correlation bins ranging from C.9 to lO'.O. We checked with an 
analysis of three and five bins that our results do not depend sig- 
nificantly on this choice. This step is followed by a minimum 
chi-square analysis and the likelihood contours in the contam- 
ination fraction parameter space are presented in Fig. Ql] The 
degeneracy between the two contamination fractions is clearly 
evident, and lower and upper limits can be estimated. 



For the lowest redshift cuts there is a strong degeneracy, 
with 0.0 < fiz < 0.6, and /1 ~ 0.01, hence a potentially large 
fraction of low redshift galaxies are expected to be at higher 
redshift. Likewise for the highest redshift cuts a potentially 
large fraction of high redshift galaxies are expected to be at 
lower redshifts. It is important to reiterate that the fraction /12 
is the number of galaxies that move from bin 1 to bin 2 as a frac- 
tion of the true number of galaxies in bin 1 . Hence, the observed 
large degeneracies in the low and high redshift cuts are a con- 
sequence of the very different occupation numbers in the two 
bins. For the intermediate redshift cut range 0.4 < z cut < 0.9 
the contamination factors are around 10% - 20%. 

The spectroscopic data in Wl allow us to directly calculate 
contamination fractions /12 and /21 for the field Wlp2p3 and 
to check whether our estimates obtained via correlation func- 
tions are reasonable. The VVDS in this field has the same depth 
as the limiting magnitude of our correlation analysis sample 
(i' AB < 24.0). We subdivide galaxies according to spectro- 
scopic (true) redshifts and determine directly contamination 
fractions with the photometric redshifts. Error contributions on 
this quantity are the Poisson noise and the redshift sample vari- 
ance in a field of 1 sq. d egree. The latter was estimated in 
van Waerbeke et al. (2006) to be 15 times the Poisson contri- 
bution. We plot our results in Fig.Qj] We see that the contami- 
nation fractions determined with spectroscopic redshifts in one 
field are very comparable to the correlation function estimates 
for the whole CARS surveys. This directly shows the validity of 
our analysis. 

5. Available data products 

We make available on request our multi-colour catalogues in- 
cluding photo-z estimates of 37 CARS fields (corresponding to 
~ 30deg 2 effective survey area after image masking). The data 
package includes object catalogues, the derived image masks 
and JPEG images to inspect colour data and extracted sources. 
The catalogues are available as FITS binary tables and a subset 
of the most important entries as UNIX-ASCII text files. 

The FITS catalogue version includes most of the original 
SExtractor keywords and for the ir meaning we refer to the 
SExtractor manual dBertinl 120031) . All these basic keys are 
measured in a SExtractor run in dual-image mode where we 
use the unconvolved /-band image for detection as well as for 
photometric measurements (see Sect. 13.3b . 

Additional keywords created in subsequent SExtractor 
runs with the PSF-matched images in the five bands for pho- 
tometric colour measurements are indicated by an additional 
_x where x € [u*,g',r',i',z']- In particular, these are the dif- 
ferent kinds of fluxes, magnitudes and corresponding error 
estimates (e.g. FLUX_MJTO_x, FLUXERR_AUT0_x, MAG_IS0.x, 
MAGERR_IS0_x, etc.); note that magnitude error estimates in 
the catalogues do not take into account systematic zeropoint 
offsets but only statistical errors due to photon noise. We es- 
timate 24 different aperture fluxes and magnitudes with di- 
ameters ranging from 4 to 55 pixels (=07744 to l(y.'23). We 
add the lcr limiting magnitudes MAG_LIM_x as described in 
Sect. 13.21 All magnitudes are provided in MegaPrime instru- 
mental AB magnitudes. We note that we did not apply any 
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0.0<z<0.3 0.3<z<0.3 0.3<z<04 0.4<z<0.5 0.5<z<0.6 0.6<z<0.7 0.7<z<0.8 0.8<z<0.9 0.9<z<1.0 1.0<z<1.5 1.5<z<20 2.0<z<4.0 




0.1 1 10 

0|_arcminj 



Fig. 10. The angular cross-correlation of objects in different photo-z slices (solid lines represent objects with i' AB < 24 and 
dashed lines represent objects with 23 < i' AB < 24): The matrix plot represents the complete CARS area and each panel shows the 
qualitative behaviour of auto- (diagonal) and crosscorrelation (off-diagonal) measurement for different redshift slices. The rows 
represent bins from z = to z = 3.9 (steps z = 0.0; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 1.0; 1.5; 2.0; 3.9). The lower left panel is 
a zoom-in to the corresponding panel of the matrix plot. Note that in the cross-correlation analysis we include all objects with 
ODDS > 0.8. See text for further details. 



magnitude correction to the catalogue entries also if our tests 
performed in Sect. lA.4l might justify them. This especially ap- 
plies for discrepancies present in the M*-band calibration of the 
W3 and W4 pointings (see Sect. IA.4l i. To allow an easy identifi- 
cation of objects with problematic photometry we add the flags 
NBPZ_G00DFILT indicating the number of filters with reliable 



photometry, NBPZJ3ADFILT giving the number of filters with 
MAGERR_IS0_x > 1.0 and NBPZ.LIMFILT listing the number of 
filters with MAG_IS0_x fainter than our formal magnitude limit 
(see Sect. [2j. Which of these three properties applies to which 
filters is encoded in additional keys. 
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Fig. 11. The estimated contamination fractions, /12 and /21 resulting from a strict pairwise cross-correlation analysis. The two 
redshift bins are defined to be 0.0 < z\ < z C ut and z cut < Z2 < 4.0. Contours indicate the 1, 2 and 3-sigma confidence regions, 
having progressively darker shades of grey. Legends in the figure give the observed number of galaxies in each bin. The data 
points with error bars denote the measured contamination fraction found for those galaxies in the VVDS Wl spectroscopic sample 
(see text for details). Note that for z cut = 1.5 and z cut = 2.0 an extended vertical scale is used in order to show the measured 
contamination. 



Furthermore, we provide a global mask key MASK which 
is for objects that do not lie inside one of our object masks 
and 1 otherwise. This key takes into account all masks from 
our object density, bright star and asteroid track analyses as 
described in Sect. 13.41 



Finally, the catalogues contain photo-z relevant quantities 
from the output of BPZ. Besides the Bayesian redshift esti- 
mate, Z_B, we include the ODDS probability, the SED corre- 
sponding to the Bayesian redshift (T_B), the corresponding^ 2 , 
the 95% confidence interval (Z_BJHN and Z_B_MAX) as well as 
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the maximum-likelihood redshift and type estimate (Z_ML and 
T_ML), which are put out by BPZ before the prior is applied. 

The most important catalogue entries are summarised in 
TableS 

6. Summary and conclusions 

We have presented high-quality five-band multi-colour data 
from 37 sq. deg. of the CARS survey. We gave a detailed de- 
scription of our data-handling procedures ranging from data 
selection to the final catalogues including a first set of pho- 
tometric redshift estimates. Our algorithms provide an accu- 
rate astrometric alignment on the sub-pixel level to extract pre- 
cise object colour information. For the large majority of our 
data the Elixir photometric information allows us to derive 
an unbiased absolute photometric calibration with a scatter of 
<x ss 0.02-0.05 on a pointing basis for g'r'i'z'; tests against the 
official TERAPIX TSSS3 CFHTLS-data release show very sig- 
nificant zeropoint offsets for four out of 93 common fields. In 
u* direct comparisons with SDSS suggest that our zeropoints 
are systematically about 0.1 mag too faint. 

We showed that our colour catalogues allow, with the help 
of spectroscopic information, the estimation o f reliable photo- 
metric redshift estimates with the method of Benitez (2000). 
In our 37 sq. deg. survey (about 30 sq. deg. in unmasked ar- 
eas) we detect about 3.9 million objects classified as galax- 
ies (SExtractor CLASS.STAR < 0.95). From those about 1.45 
million (10 - 15 galaxies per sq. arcmin) have a formally reli- 
able photo-z estimate with ODDS > 0.9 (completeness 37.2%). 
Comparing our photo-z estimates with external spectroscopic 
data we find an overall performance of cr(Az/(l + z)) ~ 0.04 - 
0.05 up to ;'ab ~ 24 with an outlier rate of 77 « 1% - 3%. We 
applied a cross-correlation analysis to qualitatively investigate 
redshift slice contamination between samples in different red- 
shift bins. It indicates significant contamination of neighbour- 
ing redshift slices with a width of Az ~ 0.1 and a dying cor- 
relation signal for bins more than Az ~ 0.3 apart. Catastrophic 
outliers occur between low-z bins and galaxies with an estimate 
of Zphot > 1.5. We performed a more quantitative analysis only 
for the case when our whole redshift sample is divided in ex- 
actly two redshift bins. With the help of spectroscopic redshifts 
from the deep part of the VVDS it reconfirms the homogeneity 
of our photo-z sample over the entire CARS area. A more com- 
plete, in-depth analysis with the correlation function technique 
will be presented in Benjamin et al., (in prep.). 

We note that the catalogues and the creation of photo-zs 
was optimised for studies in the regime < z < 1.4 and objects 
with a larger estimate should be filtered. The current catalogues 
are not suited for studies of the high-z regime such as M*-band 
drop-out searches. While the photo-z performance according to 
formal parameters is very good our estimates show a systematic 
tilt for < z < 1 (higher redshift ranges cannot be verified due 
to the lack of spectroscopic information). Our estimates are too 
high by Az * 0.1 for low z and the bias decreases linearly to 
reach about Az ~ -0.1 for z ~ I - The zero-crossing of the tilt is 
at z ~ 0.5. The mean bias is about (Az) ~ 0.03 for z < 0.5 and 
about (Az) ~ -0.03 for 0.5 < z < 1. Improved and bias-free 
BPZ photo-z estimates will be presented in Hildebrandt et al. (in 



pr eparation). Addition ally, photo-z estimates with the method 
of iBender et al.1 (1200 11) will be analysed and compared to our 
current work in Brimioulle et al. (in preparation). 

The presented catalogues mark the first step for the primary 
science goal of CARS in the CFHTLS-Wide area: The assem- 
bling of a galaxy cluster sample from low to high redshift and 
its subsequent exploitation for cosmological studies. For the 
second step in this effort, our multi-colour data are currently 
being used on several cluster detection algorithms : The Voronoi 



tessellation technique from Ra mella et al. d200ll) . the Postman 



matched filter algorithm (see Postman et al. 19961) and t he Red- 
Cluster Sequence technique (see ldadders & Yea2000h . 

To trigger a larger variety of follow-up studies we make 
available our catalogues on request. 
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Table 4. Description of the most important FITS keys in the CARS multi-colour catalogues. The ASCII catalogue version contains 
one aperture magnitude at a diameter of 1'.'86. The FITS version lists 24 aperture magnitudes for diameters from 0'.'744 to 1(X'23. 



key name 


description 


measured on 


ASCII catalogue 


SeqNr 


Running object number 




y 


ALPHA_J200® 


Right ascension 


unconvolved /'-band image 


y 


DELTA_J200® 


Declination 


unconvolved /'-band image 


y 


Xpos 


x pixel position 


unconvolved /'-band image 


y 


Ypos 


v niYpl nosition 


nnponvol vpH /'-hand imaop 

L111L-LJ11V xJ 1 V L.LI I U CLllL-i 1111 till L^ 


y 


MAG AUTO 


total /'-hand maoTntiirlp 

IV.' 1 CI J t L'tlt 1L.I 11 ILL £1111 L LILIL, 


nnponvolvpd /'-band lmatrp 

U11L^L_J11VV_/1VL^L1 / L'tlllLI 11 llclw L. 


y 


MAGFRR AUTO 


total /'-hand Tnatrmtiirlp prror 

LULCll t L'tlt 1L.I 11 ILL £1111 L LI LIC L I 1 V / 1 


nnponvol vpH /'-hand imaop 

LIllL^LJllVV-'lVL^Ll 1 U CLllL-i llllullL^ 


y 


MACISCLx" 


isophotal magnitude in x-band 


PSF-equalised x-band image 


y 


MAGERR_ISO_x 


isophotal magnitude error in x-band 


PSF-equalised x-band image 


y 


MAG_APER_x 


aperture magnitude vector in x-band 


PSF-equalised x-band image 


y 


MAGERR_APER_x 


aperture magnitude error vector in x-band 


PSF-equalised x-band image 


y 


FWHM_WORLD 


FWHM assuming a Gaussian core 


unconvolved /'-band image 


y 


FLUX_RADIUS 


half-light-radius 


unconvolved /'-band image 


y 


A.WORLD 


profile RMS along major axis 


unconvolved /'-band image 


y 


B.WORLD 


profile RMS along minor axis 


unconvolved /'-band image 


y 


THETA_J2O0O 


position angle 


unconvolved /'-band image 


y 


CLASS.STAR 


star- galaxy classifier 


unconvolved /'-band image 


y 


Flag 


SExT.ra.ctor extraction flags 


unconvolved /'-band image 


y 


FLUX_ISO_x 


isophotal flux in x-band 


PSF-equalised x-band image 




FLUXERR.ISO.X 


isophotal flux error in x-band 


PSF-equalised x-band image 




FLUX_APER_x 


aperture flux vector in x-band 


PSF-equalised x-band image 




FLUXERR_APER_x 


aperture flux error vector in x-band 


PSF-equalised x-band image 




MAG_LIM_x 


limiting magnitude in x-band 


unconvolved x-band image 


y 


Z_B 


Bayesian photo-z estimate 




y 


Z.B.MIN 


lower bound of the 95% confidence interval 


_ 




Z.B.MAX 


upper bound of the 95% confidence interval 






T.B 


best-fit spectral type^ 




y 


ODDS 


empirical odds^ 




y 


NBPZ_GOODFJ.LT 


filters with reliable photometry 




y 


NBPZ.BADFILT 


filters with MAGERR_I SO > 1.0 




y 


NBPZ.LIHFILT 


filters with MAG_IS0_x > MAG_LIM_x 




y 


MASK 


global mask key^ 




y 



8 x 6 [u, g, r, i, z] 

b Ell= 1, Sbc= 2, Scd= 3, Im= 4, SB3= 5, SB2= 6, plus two interpolated types in colour-redshift space between each pair of these basis 
templates. Intermediate best-fit templates are indicated by a floating point number for T_B. 
c integrated probability inside an interval which is such that it contains 95% probability for a single Gaussian 
d unification of the different masks described in Sect. [34] 
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Appendix A: Details on the data handling of CARS 
observations 

In this appendix we give a more detailed description on our data 
handling procedures of the CARS survey. The excellent Elixir 
preprocessing of the CARS data and the available meta-data in- 
formation (see below) allowed us to build up a complete sur- 
vey pipeline starting from data retrieval up to final co-added 
science images. The construction of our data-processing and 
many choices for our data-handling as described below were 
driven by the following two requirements: 

1 . The system allows a 100% automatic processing of the data 
with the need for manual intervention only at the final veri- 
fication stage of co-added science images. This particularly 
forbids manual passes through all individual Elixir im- 
ages either to visually grade data or to remove/mask arte- 
facts. We needed to automatically reject problematic ex- 
posures or parts of them from the whole analysis or we 
need to deal with remaining defects at the level of the final 
science images. For instance, our automatic satellite track 
removal module reliably detects and removes about 95% 
of all bright satellite trails. If necessary, the remaining 5% 
need to be masked manually in the final science images. 

2. We want to independently and incrementally process indi- 
vidual pointings as soon as the full five band coverage of a 
particular area is becoming publicly available. 

We completely achieved the first goal and the complete data 
processing of the CARS data is done by one of the authors with 
two computers (a double processor Athlon2800+ with 4 GB 
Virtual memory and a quad processor/dual core AMD Opteron 
885 with 1 1 GB of RAM) and a total disk storage capacity of 
about 10 Terabytes. The second goal could not be met for the 
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photometric calibration of several fields and we needed to use 
information from adjacent pointings to obtain an absolute flux 
calibration (see below). 

Most of our algorithms to process optical dat a from 
multi -chip cameras were d escribed in lErben et aL I d2005h 
and iHildebrandt et all d2006h in the context of GaBoDS data 
(see e. g. ISchirmer et alJ 120031) from the 8-chip instrument 
WFI@MPG/ES02.2m. We therefore limit the discussion to pe- 
culiarities of the CARS data, necessary pipeline upgrades due to 
the four times larger field-of-view of MegaPrime and quality 
assessments of our final science images. 

A .1. Data preselection and retrieval 

As described in Sect. [2] the starting point of the current 
CARS data set are the Elixir preprocessed images from 
the CFHTLS-Wide Survey. Besides the images, compre- 
hensive information on the current status of the Survey 
and the observed data is available in the form of a CFHTLS 
exposure catalogue (see |http : //www . c fht . hawa iTTj 
edu/Science/CFHTLS-DATA/exposurescatalogs Thtml] ). 
From the CFHTLS-Wide we preselect all pointings which are 
publicly available at 18/01/2008 and which have observations 
in the complete u*g'r'i'z' filter set. 

For each survey-image a Service Observer quality 
flag ranging from 1-5 is available. A '1' states that the ex- 
posure was obtained within survey specifications and has no 
obvious defects. A '2' means that one of the predefined spec- 
ifications for that exposure (seeing, sky transparency or moon 
phase) was out of bounds. A flag of three or higher indicates 
poor observing conditions or other severe defects such as track- 
ing problems during the exposure. Only images with flags ' 1 ' 
or '2' enter our processing. We visually inspected in total 300 
Elixir preprocessed MegaPrime exposures (60 in each filter) 
to verify the suitability of this quality assessment for a blind 
and automatic preselection of CARS data. On the other hand we 
did not check whether a subset of the images with higher flags 
still could be included in our survey. 

We use the aforementioned information and the possibility 
to request CADC files and data products directly within pro- 
grams or shell script^] to automatically retrieve the images 
of interest. For the current work we transfered in total 1246 
MegaPrime CFHTLS-Wide images from CADC. 

A. 2. Data preprocessing and weight image creation 

The Elixir preprocessing includes all necessary opera- 
tions to remove the instrumental signature from raw data. 
The data on which we start our analysis are bias-corrected 
and flat-fielded. Moreover, fringes are removed in i' and z' 
observations, permanent bad CCD pixels are marked and 
all images are corrected for photometric non-uniformities 
across the MegaPrime field-of-view; see the WWW pages 
|http : //www . cfht .hawaii . edu/Science/CFHTLS-DATA/| 
datapro cessing . html and |http : //www . cfht .hawaiiTj 



edu/Science/CFHTLS-DATA/megaprimecalibrat ion. 
html for a more detailed description of the Elixir processing. 

The visual appearance of the Elixir processed data is very 
good. Only fringe residuals are observed for parts of the z' data. 
As discussed within Fig. [2] individual chips of certain expo- 
sures might not contain useful data. The first step of our own 
processing is therefore to identify problematic chips by con- 
sidering pixel statistics and to mark them as unusable. More 
precisely, we exclude chips with the following defects from 
any further analysis: (1) The pixel value at the lower quartile 
of the chip pixelvalue-distribution is 10 or lower. This means 
that large fractions of the chip contain zeros; (2) More than 3% 
of the pixels in a chip are saturated. This means that a consid- 
erable chip-area is contaminated by a very bright star which 
would most probably lead to problems in the later astrometric 
calibration. Furthermore, at the level of science analysis such 
areas would be excluded anyway; (3) To astrometrically cali- 
brate our data we first tie the i' observations to the USN0-B1 
catalogue and then extract from the stacked V image a deeper 
astrometric reference system for the other colours. This is dis- 
cussed in more detail in Sect. IA.3I below. For this reason we 
exclude from u*g'r'z' all chips which have been identified as 
bad in the corresponding /'-band data. In the following we up- 
dated all our THELI pipeline modules to smoothly handle an 
arbitrary geometry of usable chips within a CCD-array. 

To prepare the extraction of object catalogues for astromet- 
ric calibration we create for each chip a corresponding weight 
map which contains information on bad pixels and the relative 
noise properties of the image pixels. The steps o f our weight 



image creation are described in detail in Sect. 6 of Erb en et al 



(2005). For the CARS data we updated and expanded our proce- 
dures as follows: 

1 . For the WFI@MPG/ES02.2m camera, the starting point of 
our weight maps was a normalised flat-field image. A flat- 
field maps the relative sensitivity of image pixels within 
a CCD array and allows us to take into account associated 
pixel noise variations during object extraction. For the CARS 
data we neglect this effect and start the weight creation with 
a flat image with a pixel value of ' 1' on the whole array. We 
verified that this simplification has no significant effect on 
our object catalogues of single frames and of the final co- 
added science images later-on. It allows us to store weight 
images very efficiently and with significantly reduced hard- 
disk space. 

2. Permanent bad pixels of the CCDs are marked in the 
Elixir processed science frames by a pixel value of '0'. 
This information is transfered to our weight maps. For 
WFI@MPG/ES02.2m data defect pixels had to be iden- 
tified with dark frames and/or flat-field images. 

3. We used to visually identify and to mask bright satellite 
tracks which must be excluded from the object extrac- 
tion and co-addition process. To process the CARS data 
we developed an automatic track detection and mask- 
ing tool based on Hough transfor m techniques (see e. g. 
Duda &Hartlll972l: IVandamdliool . To reliably find real 



11 see | http : / /wwwl . cade- ccda . hia- iha . nrc- enre . gc . ca/ 
getData/doc 



tracks and to reject spurious detections due to bright stars 
and extended objects we use that a satellite typically con- 



20 



T. Erben et al.: CARS - Five-band multi-colour data from 37 sq. deg. archival CFHTLS observations 



taminates several chips on the MegaPrime mosaic. When a 
candidate is found on a particular chip we check for de- 
tections on expected positions in other detectors. Tracks 
which are found in two or more chips are masked on all 
CCDs crossing their path including an additional margin 
of one CCD on both sides of the track. In this way also 
detectors on which a track cannot be detected individually 
are appropriately covered by an image mask. From man- 
ual inspection of 60 MegaPrime exposures it is found that 
our implementation correctly detects and masks more than 
95% of all bright satellite tracks. The failures can mostly be 
traced to either dashed, non continuous satellite tracks or to 
short ones at the edges of the mosaic. We observed only a 
handful false-positive detections in very special configura- 
tions such as extended and bright object chains over chip 
boundaries. 

Pixels attached to identified tracks are set to zero in the 
corresponding weight maps. 

A. 3. Astrometric calibration 

With the preprocessed Elixir images and the weight maps, 
the astrometric calibration of CARS data sets and the associated 
quality assess ment follows very closely the procedures outlined 
in Sect. 5 of lErben et all d2005l) : 



1. 



SExtractor is used on all images to extract sources with 
at least 5 pixels having 5<x above the sky -backgroun d. 



2. The Astrometrix 



McCracken et al 



programme (see 
2003) is run on the 



Radovich 2002; 



/'-band to deter- 
mine a third order astrometric solution for each individ- 
ual chip. We use t he USN0-B1 standard-star catalogue (see 
Monet et al. 2003) as our astrometric reference frame. We 



co-add the /'-band data (see below) and extract a high S /N 
object catalogue (sources with at least 20 pixels of 20<x 
above the sky-background noise) from it. This catalogue 
is used as reference for the astrometric calibration of the 
u'g'r'z' data. It allows us to use a dense reference cata- 
logue with high positional precision. This turned out to be 
essential to map robustly significant higher-order astromet- 
ric distortions of MegaPrime. For this instrument, the con- 
tributions of second- and third-order terms to the astromet- 
ric correction are up to 2'.'0. As a comparison, for the four 
times smaller WFI@MPG/ES02.2m camera these values 
are on the order of 075; see also Fig. IA.ll We use the /'-band 
data as reference because (1) the individual exposures are 
already reasonably deep and hence can be well calibrated 
with the USN0-B1 sources and (2) with seven dithered ex- 
posures the co-added /'-band typically has the best cover- 
age and filling-factor of the MegaPrime area with its chip 
gaps in individual exposures. 

For our scientific objectives, weak gravitational lensing, 
and multi-band studies with photometric redshifts it is es- 
sential to obtain a very high internal astrometric accu- 
racy in the lensing band (/' in our case). As discussed in 



Erben et a l. (2005) higher order object brightness-moments 
(everything above the zeroth order moment, i. e. the ob- 
ject flux) are significantly changed if individual frames of 



the WFI@MPG/ES02.2m camera are aligned with an ac- 
curacy > 0.5pixel (pixel scale 07238). For our CARS data 
we reach an accuracy of about 07 03-07 04, i. e. about 1 /5th 
of a MegaPrime pixel. This was tested by comparing object 
positions (with /ab < 20) from astrometrically corrected in- 
dividual frames with their cousins in the final co-added im- 
ages. Results for the field Wlmlp2 in the /' filter are shown 
in Fig. IA.2I In a similar way the inter-colour alignment of 
different filters per pointing is tested. From the final co- 
added images we compare object positions of high S /N 
sources (/ab < 20). For all pointings we reach an align- 
ment between the colours below one pixel. We note how- 
ever that the inter-colour alignment between g'r'i' is bet- 
ter (typically 0.5 pixels) than between i'u*z' (between 0.5 
and 1 pixel). This is expected because individual frames 
from u* and z' have fewer high S /N sources for astrometric 
calibration than the intermediate filter bands. We show re- 
sults from the /' - g' and /' - u* comparisons of Wlmlp2 in 
Fig.[0 



A.4. Photometric calibration 

The Elixir preprocessed images come with all neces- 
sary meta-data to translate pixel counts to instrumental AB 
magnitudes; see |http : //www . cfht . hawaii . edu/Science^ 
CFHTLS-DATA/megaprimecalibration.html#P2 for a full 
description of the Elixir procedures to derive photometric 
parameters. In addition, the CFHTLS exposure catalogue 
contains a flag whether an image was taken under photomet- 
ric sky conditions or not. With this information we try to derive 
a photometric zeropoint for each colour of each pointing (we 
call this a set in the following) w ith the two-stage process de- 
scribed in Hildebrandt et al. (2006): 



1. We use Photometrix to bring all individual images to the 
same flux scale by estimating the magnitude differences of 
overlap sources. This gives us for each image i a relative 
zeropoint ZP K \i which tells us the magnitude offset of that 
image w.r.t. the mean relative zeropoint of all images, i.e. 
we demand £;ZP re y = 0. Note that this procedure rela- 
tively calibrates data obtained under photometric and non- 
photometric conditions. An absolute flux scaling can now 
be obtained from the photometric subset. 

2. For images being observed under photometric conditions 
we calculate a corrected zeropoint ZP con j according 
to ZP mn j = ZP + am ■ ext + ZP te \j, where am is the air- 
mass during observation, ZP the instrumental AB zeropoint 
and ext the colour dependent extinction coefficient. The last 
two quantities are part of the provided Elixir meta-data. 
For images obtained under photometric conditions, the rel- 
ative zeropoints compensate for atmospheric extinction and 
the corrected zeropoints agree within measurement errors 
in the ideal case. We later use the mean (ZP corrj ) as zero- 
point for our co-added images. 

Note that we can determine an absolute zeropoint only if at 
least one exposure of a given set was obtained under photo- 
metric conditions. For the CARS data this is not the case for 
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Fig. A.l. Higher order MegaPrime distortions: The plots show 
the difference in object position after a first-order astromet- 
ric alignment, i. e. corrections for linear shifts and rotations 
and a the full third-order astrometric solution estimated by 
Astrometrix. The sticks in the upper plot indicate the posi- 
tional displacement vector between the two solutions and the 
lower plot gives the absolute displacement numbers. We note 
that the second- and third-order terms contribute very signifi- 
cantly to the solution (see text for details). 

12 sets and we estimated a zeropoint for those after image 
co-addition by flux comparison with objects of adjacent, cali- 
brated pointings. As discussed in Hildebrandt et al. ( 20061) the 
corrected zeropoints offer a good opportunity to verify 
the quality of absolute photometric calibration. If a field is ob- 
served over an extended period, the comparison of zeropoints 
estimated from different nights gives a robust indication on the 
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Fig. A.2. Internal Astrometric alignment of /'-band data from 
Wlmlp2: Shown are the differences of sources in 7 dithered i'- 
band observations after astrometric calibration with the coun- 
terparts in the co-added image. The 7 individual exposures 
were obtained with a dither pattern spanning about 4570 in 
Ra and 18070 in Dec to cover the MegaPrime chip gaps. 
The plot covers the complete MegaPrime area of about 1 
sq. degree. The thick solid lines mark the region containing 
68% of all points and are at ARa = -0.003+jJ jjg arcsec and 
ADec = -0.002+QQjg arcsec. Dashed lines show the corre- 
sponding area for 90% of all points (ARa = -0.003*9 9?? arcsec 
and ADec = -0.002+°°^) arcsec. Only 1 out of 10 points is 
shown for clarity of the plot. 



long-term stability of photometric instrument properties and on 
the calibration process itself. 

However, the CARS fields from the CFHTLS-Wide Survey 
were observed mostly in a compact period during a single 
night and hence do not allow for this test directly. To perform 
this important quality control we consider observations of the 
CFHTLS-Deep Survey. This part of the CFHTLS continously 
observes four one square degree fields with the goa l to detect 



Supernovae and to measure their light curves; see lAstier et al 



(2006) for more details on this survey. Because we noticed 
severe problems with our u* flux calibration for W3 and W4 
later-on (see Sect. lA.6l below) we investigated CFHTLS-Deep 
images which were publicly available at 01/01/2008 and 
which had the photometric flag in the CFHTLS exposure 
catalogue. We studied the long-term evolution of the cor- 
rected zeropoints from June 2003 to December 2006. Results 
for the fields Dl, D3 and D4 are shown in Fig. IA.4I For the 
Dl area we studied in detail photometric stability in all five 
filter bands. Fig. IA.4I shows that the photometric calibration 
of this field over time is very consistent with formal standard 
deviations of only up to about 0.03 mag. Note however that the 
peak differences of magnitude zeropoints span more than 0. 1 
mag! For the M*-band observations of D3 and D4 we observe a 
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Fig. A.3. Astrometric alignment of different colours of CARS 
field Wlmlp2: We show differences in object positions from 
different colours u*g'i' of the CARS pointing Wlmlp2. The 
plot shows that the different bands are well aligned with sub- 
pixel accuracy although slight trends in the residuals with u* 
are visible. Solid (dashed) lines enclose areas containing 68% 
(90%) of all points. They are at x v - % = -0.006^(+o.iO) 
arcsec, y,< - y g > — -Q.OH^q^q^) arcsec and xy - x u - - 
-0.009!^C°;|6) arcsec, y-,- y~, '= -0.03^ 2 (^ arcsec. 
See text for further details. 

considerably larger scatter with extremely low values for the 
corrected zeropoints (too high CARS magnitude zeropoints), 
from around April 2006 to November 2006. This hints to a 
calibration problem of M*-band CFHTLS data in that period 
and the data suggest a necessary correction of about -0.2 
mag for M*-band CARS observations from spring to fall 2006. 
Due to its visibility Dl has no observations in that period. 
We verified that the D4 results for the colours g'r'i'z' are 
similar to those in Dl and hence the problem seems to be 
confined to u*. The se discrepanci es in the u* calibration are 
also documented in iGwynl J2008h and on TERAPIX WWW 
pages describing the official T00(54 CFHTLS data release, 
(http ://terapix. iap. fr/article.php?id^article=713) 

The Deep data allow us also to check homogeneity 
and reproducibility of our photometric calibration over the 
MegaPrime field-of-view. In each colour we created three in- 
dependent co-added images from Dl. Each stack contains five 
images obtained under photometric conditions in November 
2003, 2004 and 2005. From the different stacks in each colour 
we match bright (17 < m < 20) sources whose positions agree 
to ff.'S or better. Magnitudes are compared and mean offsets 
and standard deviations are estimated. Table I A. 1 1 summarises 
the results. Except for the z'-band the total magnitudes agree to 
better than 0.04 for these stacks. The scatter around these abso- 
lute offsets is around 0.02 mag for u*g'r'i' and about 0.03 -0.04 
mag for z'\ we attribute the higher value in z' to fringe residuals 



Fig.A.4. Long-term studies of corrected zeropoints in 
CFHTLS-Deep data: We analysed all publicly available 
photometric data from the CFHTLS-Deep fields Dl (Ra: 
02:26:00; Dec: -04:30:00; all colours), D3 (Ra: 14:19:28; 
Dec: +52:40:41; M*-band) and D4 (Ra: 22:15:31; Dec: 
-17:44:06; «*-band). The panels show the long-term evo- 
lution of corrected photometric zeropoints in these 
fields from 25/06/2003 - 29/12/2006. On the measured zero- 
point distributions we performed an iterative 3cr clipping to 
exclude obvious outliers. The quoted GaBoDSID is a running 
number counting the nights from 3 1/12/1998. The vertical lines 
in the D3 and D4 panels show the period of W3 and W4 u* obser- 
vations for which we observe larger discrepancies in compar- 
isons with SDSS magnitudes; see text for further details. 

in this band. These values give us an estimate on the internal 
photometric accuracy of our data, i. e. on the error propaga- 
tion of inaccuracies in our photometric calibration procedures 
to final magnitude estimates. We consider them as upper limits 
because possible errors in the determination of the photometric 
superflat within Elixir contribute to the quoted numbers. We 
note that we obtain consistent errors when comparing the Deep 
stacks with the CARS pointing Wlp2p3. Their distance on the 
sky is 5'.0 in Ra, 18'.0 in Dec and we can compare fluxes from 
objects which fall in different areas of the MegaPrime mosaic. 
The test between the Deep stack from 2003 and Wlp2p3 yields: 
Am* =0.003 +0.022; Ag' = 0.014+0.014; Ar' =0.053+0.012; 
A? = -0.019 + 0.014; Az' = -0.034 + 0.034. 

A.5. Image co-addition 

After the photometric calibration we check whether the indi- 
vidual exposures of a given set were obtained under varying 
photometric conditions which is indicated by a large range of 
relative zeropoints. Low-value outliers point to images which 
were observed under unfavourable sky-conditions w.r.t. the rest 
of the set. We estimate the median med(ZP K ij) of the relative 
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Table A.l. Comparison of magnitudes from three independent 
stacks of CFHTLS-Deep Dl data in each colour: The upper 
row shows magnitude differences between stacks of 2003 and 
2004, the lower one between data from 2003 and 2005. See text 
for details. 



W1p2p2 (i'-bond) 



Au' x 100 


Ag' x 100 


Ar' x 100 


Ai' x 100 


Az' x 100 


-3.3 ± 1.8 


1.3 + 1.5 


-0.5 ± 1.3 


-1.4 ±2.0 


4.6 ±3.1 


-2.6 ±2.0 


3.3 ± 1.6 


0.6+1.6 


-0.0 ± 2.0 


-7.3 ±3.8 



zeropoint distribution and sort out exposures with a relative ze- 
ropoint of med(ZP le \j) - 0.1 or smaller. Affected are the five 
CARS sets Wlplp2-i ' , Wlp3m<5-i ' , Wlp4m0-i ' , W4mlml-z ' 
and W4mlm2-z ' . At this stage we also reject short calibration 
exposures (f exp < 100s) from further processing. 

Finally the exposures belonging to a set are sky- 
subtracted with SExtractor and co-added with the Swarp 
programme (see iBertinl 120081) . We use the LANCZOS3 ker- 
nel to remap original image pixels according to our astro- 
metric solutions. The subsequent co-addition is done with 
a statistically optimal weighted mean which takes into ac- 
count sky-background noise, weight maps and the relative 



photom etric zeropoints as described in Sect. 7 of Erb en et al 



(2005 ). As sky projection we use the TAN projection (see 
Greisen & Calabrettal Eool and all colours from a specific 
pointing are mapped on the same pixel grid. The origins of 
the TAN projection for each pointing are those defined for 
the CFHTLS -Wide survey (see http : //terapix . iap . fr/| 
cplt/oldSite/Descart/summarycfhtlswide .html). An 
example for a final co-added image is shown in Fig. [2] 

After co-addition we cut all images to a common size of 
21k x 21k which covers areas with useful data for all CARS 
pointings. We visually check each product for obvious defects. 
Because no pixel rejection takes place in our weighted mean 
stacking the co-added images show some remaining artefacts. 
We observe faint satellite tracks which were not detected and 
removed by our automatic satellite track masking tool (see 
Sect. \A.2i , short asteroid trails and warm pixels. Within the 
subsequent catalogue creation process we tried to mask these 
defects on the basis of the /'-band of each pointing (see Sect.|3]l. 

A first rough check for the photometric calibration of each 
co-added image is done with the help of galaxy number counts. 
In all colours we estimated reference counts from the CFHTLS 
Dl field and compare those to our CARS Wide data. The ref- 
erence counts are about two magnitudes deeper than those 
from our co-added images. Each of the one square degree 
fields yields a robust estimation of these counts and allows 
us to quickly spot photometric calibrations with obvious prob- 
lems (Am 0.2). Galaxies are selected by the SExtractor 
CLASS_STAR parameter (CLASS.STAR < 0.95); see Fig. |A31 
for an example of the field Wlp2p2. With multi-colour obser- 
vation for all fields we also can c ompare colours from stellar 
sources with predictions from the IPicklesI ( 1998 ) library. We 
select bright, unsaturated stars by magnitude (17 < i' hB < 22) 
and by CLASS_STAR > 0.95. With data in u*g'r'i'z' we plot ten 
possibl e colour-colour combinations against the Pickles stellar 
library dPickleslll998l) which allows us, similar to the galaxy 



^ 3 



— , 1 1 , , 1 , 1 , 1 , , r- 

goloxy no.= 87627 

linefit for mag in [17.B;23.3] with slope = 0.360 

crosses: W1p2p2 galaxy numbercounts 
diamonds: D1 comparison numbercounts 
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Fig. A.5. /'-band numbercounts for the field Wlp2p2: galaxies 
are selected with CLASS_STAR < 0.95 to allow for a first crude 
check of our magnitude zeropoints with galaxy-number counts. 

number counts, to identify grossly inaccurate zeropoints with 
Am 0.1 - 0.2. See Fig. lA.6l for an example track of Wlp2p2. 

A much more rigorous and accurate test for the photometric 
quality of our data is given by direct comparison with external 
and well calibrated data sets. This is the topic of the following 
sections. 

A.6. Comparison of CARS data with SDSS photometry 

The overlap of all three CARS patches with the SDSS (see 
Sect. I3 allows us a direct comparison of object fluxes with 
Sloan photometry. To convert the instrumental MegaPrime AB 
magnitudes from stellar objects to the SDSS system we use the 
following transformation formulae: 



(A.l) 



The relations for g'r'i'z' were determined within the CFHTLS- 



M AB 


= WSDSS 


- 0.241 


• ("SDSS 


- gSDSs) 


Sab 


= gSDSS 


-0.153 


• (gSDSS 


- ?"sdss) 


r AB 


= '"SDSS 


- 0.024 


(gSDSS 


- 7"SDSS) 


'ab 


= 'SDSS 


- 0.085 • 


(>"SDSS " 


" 'SDSS) 


i 

Z AB 


= ZSDSS 


+ 0.074 


('SDSS " 


" ZSDSS) 



Deep Supernova project (see http : //www . ast ro . uvic ■ j 
|ca/~pritchet/SN/Calib/ColourTerms-2006Junl9/^ 
index .html#Sec®4); the u* transformation comes 
from the CFHT instrument page (see |http://cfht. 
[hawaii . edu/I nstruments/Imaging/Meg aPrime/| 
generalinf ormation .html]). 

For all the following photometric comparison studies we 
extracted single frame photometric catalogues from all CARS 
images and we use the SExtractor MAG_AUT0 estimate 
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W1p2p2 




Fig.A.6. u* - r' vs. r' - V colour-colour diagram in the 
field Wlp2m2. The dots represent the measured colours of stars 
(SExtractor CLASS.STAR > 0.95 and / AB < 20; N = 1089) 
and the star s ymbols are colours of stars from the library of 
Pickles! dl998l) . 



throughout, i. e. here we do not use the multi-colour cata- 
logues described in Sect. [3] We compare magnitude estimates 
from sources classified as stars in the SDSS and having a CARS 
MAGJUJTO estimate of 17 < m < 20. Representative results 
of our SDSS comparisons are shown in Fig. IA.7I A complete 
listing of the measured magnitude offsets and dispersions can 
be found in Table \A2\ We note a stable calibration in g'r'i'. 
For nearly all fields the mean offset in these filters is well be- 
low 0.05 mag and the transformation relations from eq. (IA. 1 b 
are valid with o" ? , « 0.02 - 0.04 in the magnitude range 
17 < m < 20. For z' the mean offset reaches up to 0.08 
mag and also the dispersion broadens to ay « 0.04 - 0.07. 
Larger disagreements are observed for the u* filter. For the 
three Wl fields with SDSS overlap we measure a consistent off- 
set of u' - msdss + 0.241 • (msdss - gSDSs) ~ 0.1. Two of the 
u* fields (Wlp3m0 and Wlp4mQ) were observed on 01/01/2006 
and 02/01/2006 and the third one, Wlplml, on 13/12/2006 and 
hence we obtain this result for different calibration periods. Our 
long-term zeropoint analysis of Dl for which «*-band data have 
been obtained during December 2006 does not indicate larger 
systematic calibration offsets. Hence, at the current stage, we 
have no explanation for this high, consistent offset between 
our Wl u* fluxes and the SDSS magnitudes. Even consider- 
ably larger absolute offsets (up to 0.3 mag) are observed for 
all CARS u* pointings of W3 and W4. However, as was discussed 
in Sect. I A. 41 systematic zeropoint offsets for the u* calibration 
are observed from April 2006 to November 2006 and all W3 
and W4 u* band observations were obtained just in that pe- 
riod. The data presented in Fig. lA.4l suggest a necessary correc- 
tion for the M*-band zeropoint of about 0.2 mag which would 



make the observed offsets consistent with the Wl results. For 
m*, the dispersion with the transformation in eq. dA.U is about 
ov * 0.03 -0.06 (17 <u* <20). 
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Fig.A.7. Comparison of CARS fluxes and SDSS magnitudes: 
We show magnitude offsets of CARS data with SDSS overlap. 
The upper row shows the comparison for all five bands of the 
field Wlplml. The dashed line in the plots indicates the zero 
offset and the solid line the observed mean difference; see text 
for a discussion of the results. 



A.7. Photometric comparison of CARS images with 
publicly available CFHTLS data 

We further check the quality of our photometric calibration 
by comparison with publicly available CFHTLS data. We 
use the TERAPIX T00O3 data release^ (see | http://wwwl. 
cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadcbin/cfht/ 
wdbi . cgi/cfht/tpx_f ields/f ori^ and |http : //terapix .~| 
iap . fr/rubrique .php?id-rubrique=2 88) which consists 
of all available CFHTLS-Deep and CFHTLS -Wide observa- 
tions until 12/09/2005. It overlaps with our data on 24 sq. 
degrees with the full five-colour coverage on nine sq. degrees. 
Furthermore, Stephen Gwyn considers public CFHTLS 
data within his MegaPipe project which aims at providing 
calibrated and co-adde d data from the complete MegaPrime 
archive at CADC (see iGwvnl 120081) (CARS overlap on 22 sq. 
degrees with five-colour coverage on two sq. degrees). All 
three data sets start their processing from the Elixir images 
but each pipeline derives the photometric calibration with 
different software modules and by including different internal 



12 At the time of writing the most recent TERAPIX T0804 release 
(with a 35 sq. degree five-colour CARS overlap) was not publicly avail- 
able. 
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and external data sets. Hence this comparison gives us another 
check on the accuracy and limitations of our algorithms: 

- The T0003 data are processed on a patch basis, i. e. to de- 
rive a photometric solution all available information from a 
CFHTLS-Wide patch are considered simultaneously. The 
global photometric analysis takes into account overlap 
sources from adjacent pointings and also allows modest 
variations of the derived Elixir zeropoints to better ensure 
a consistent solution on the complete patch. In contrast, the 
CARS data are treated on a strict pointing-by-pointing basis. 
Other differences in the processing which might influence 
direct flux comparisons between T0003 and CARS are: 

1. To select suitable Elixir images for further process- 
ing TERAPIX does not rely solely on the quality flag in 
the CFHTLS exposure catalogue but each image is re- 
graded. Hence, the composition of image stacks might 
be different for some pointing/colour combinations. We 
did not investigate this in detail. 

2. The T0003 stacks are created with a median co-addition 
whereas CARS uses a weighted mean statistics. Satellite 
tracks in individual frames are not masked before co- 
addition in the T00Q3 processing. 

- The MegaPipe project directly uses the SDSS to photomet- 
rically calibrate MegaPrime data on a pointing basis. For 
observations which overlap with Sloan, the relations from 
eq. ( IA.U are used to derive zeropoints for the images. For 
observations outside the Sloan area, the procedure is sim- 
ilar to ours. If data are obtained under photometric condi- 
tions, the Elixir calibration is used. Otherwise, a calibra- 
tion with adjacent pointings having photometric informa- 
tion is tried. 

Similar to the T0003 processing all MegaPipe images are 
rechecked manually for their suitability to be processed 
further and the final stacking is done with a median co- 
addition. 

The results of our flux comparisons with 93 TO0O3 and 62 
MegaPipe fields are detailed in Table lA.2l 

We note in general a very good agreement between our 
calibration and that from TERAPIX T0003. For all but four 
pointings the discrepancy is less than 0.04 mag. Notable dif- 
ferences occur for the stacks Wlp2p3-r' (wcars - wtooo3 = 
-0.063 mag), Wlmlp3-M* (wcars - '«tooo3 = -0.18 mag), 
Wlmlp3-z' (mcARS - '«tooo3 = -0.12 mag) and Wlmlp2-w* 
(twcars - '«tooo3 = -0.12 mag). As discussed in Sect. IA.4l the 
field Wlp2p3 overlaps with CFHTLS-Deep 1 and we can con- 
firm an offset of about 0.05 mag between the CARS Wlp2p3-r' 
stack and corresponding Deep data. All four individual images 
contributing to the CARS image have been obtained under pho- 
tometric conditions which is confirmed by a very narrow dis- 
tribution (about 0.01 mag) of relative zeropoints. At the current 
stage we do not have a conclusive explanation for the observed 
discrepancy in this field. We note that Wlp4p3-r' for which 
we observe no discrepancy (wcars - '«tooo3 = -0.007) has 
been observed in the same night (23/08/2003) as Wlp2p3-r'. 
Furthermore, these two sets share the same photometric cali- 
bration data. 



For the other three cases with a fairly large magnitude shift 
of more than 0. 1 mag all the science frames were obtained un- 
der non-photometric conditions with a flux absorption of about 
0.2 mag which probably leads to larger errors in the estimation 
of fluxes and relative zeropoints. The images were absolutely 
calibrated with one short exposed image obtained under photo- 
metric conditions. 

The direct comparison of CARS and the MegaPipe images 
shows considerably larger scatters. We investigated in more de- 
tail the case of MegaPipe W3mlm2-r' which shows a magni- 
tude offset of nearly 0. 1 mag w.r.t. CARS, T0Q03 and the SDSS. 
It turned out that an image obtained under unfavourable pho- 
tometric conditions was included in the calibration and stack- 
ing process although it should have been rejected. The median- 
stacking of heterogeneous data (MegaPipe does, by default, 
not reject very short calibration exposures as CARS and T8003) 
and problematic images that slipped through the grading pro- 
cess probably account for the observed scatters in other cases 
(S. Gwyn, private communication). 

A.8. Photometric accuracy of CARS data - Summary 

In the preceding sections we evaluated the internal and external 
photometric quality of our data. The results can be summarised 
as follows: 

- We evaluated the internal photometric accuracy of our co- 
added data with observations from the CFHTLS-Deep sur- 
vey. From Dl we constructed in each colour three inde- 
pendent stacks which contain data from 2003, 2004 and 
2005 and compared fluxes from overlap sources. The CARS 
set Wlp2p3 which was obtained in 2003 and is offset to 
Dl is included in these tests. Around some absolute off- 
sets, the magnitude comparisons show internal scatters with 
o~int,u'g'r'i' ~ 0.01 - 0.02 in u'g'r'i' and about o- intz i » 
0.03 - 0.04 in z' and we quote these values as internal mag- 
nitude uncertainties over the MegaPrime field-of-view. 

- The accuracy of the absolute photometric calibration is 
primarily tested with a comparison to the Sloan Digital 
Sky Survey. The available Elixir pre-calibration allows 
us to obtain an absolute photometric accuracy of about 
Vabs,g'r'i' ~ 0.01 - 0.04 mag in the g'r'i' bands. Unbiased 
results are also obtained for z' with an accuracy of <r a h s ,-< ~ 
0.03-0.05 mag. At the current stage we obtain a systematic 
bias of about 0.1 mag for the M*-band. This holds directly 
for our data in Wl. For W3 and W4 we arrive at the same 
conclusion if we take into account systematics revealed by 
our zeropoint study of CFHTLS-Deep data. Given this re- 
sult we quote the zeropoint uncertainty in the M*-band with 
(r a bs,u> ~ 0.15 mag. 

- Because we process our data on a pointing basis we also 
need to calibrate our images set by set relying on Elixir 
meta-data only. We generally do not take into account in- 
formation from adjacent pointings. TERAPIX T0003 data 
are treated with a more sophisticated procedure using all 
available information to simultaneously calibrate data on a 
patch-wide basis. Our direct comparison shows that both 
pipelines lead to very comparable results with a small av- 
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erage magnitude scatter of about 0.02 mag. However, we 
observe four significant outliers (out of 93 common CARS- 
T0003 sets) with magnitude offsets of 0.05 - 0.18 mag. 
Unfortunately no other external comparison is available for 
these fields. If we take the conservative approach to at- 
tribute these offsets to inaccuracies in our calibration and 
if the current CARS data set is representative less than 5% 
of our images severly suffer from a non-optimal photomet- 
ric calibration procedure. 



A.9. Detailed CARS data quality information 

In Table lA.2l we provide detailed information on the character- 
istics of each CARS data set. It contains the effective area of each 
field after image masking (see Sect. 13.4k the number of individ- 
ual images contributing to each stack, the total exposure time, 
the limiting magnitude as defined in Sect. [2] magnitude com- 
parisons with Sloan, the TERAPIX T0003 and the MegaPipe 
releases as described in Sect. I A. 61 and Sect. IA.7l the measured 
image seeing and special comments. The comments field lists 
notable defects originating from the data itself or from our re- 
duction process. We do not list defects of astronomical origin 
(e.g. very bright stars, external reflections) or problems which 
are present in a large number of images (e.g. faint satellite 
tracks, asteroid tracks, residual warm pixels, low-level fringe 
residuals which are visible in most of the z' images). We use 
the following abbreviations: 

- no ch. XX: The stack contains no data around chip posi- 
tion^) XX. We number the MegaPrime mosaic chip from 
left to right and from bottom to top. The lower left (east- 
south) chip has number 1, the lower right (west-south) chip 
number 9 and the upper-right (west-north) chip number 36. 
Note that this labeling scheme differs from that used at 
CFHT 

- fr. res.: The co-added image shows significant fringe resid- 
uals. 

- m. ZP: The zeropoint for this image was obtained manually 
by comparing object fluxes from the image with adjacent, 
photometrically calibrated pointings; see Sect. IA.4l 

- sat. tr.: The co-added image shows a bright satellite track 
which was not masked by our track detection module. 
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Table A.2. CARS data quality overview: The first column lists: CARS field-naming convention; CFHTLS field-naming convention; effective field 
area after image masking. The seventh column contains magnitude comparisons of CARS fields with those from TERAPIX T0003 (indicated 
by (T)) and MegaPipe (indicated by (M)). Magnitude offsets are always given as m C ARS - m other . See the text for more details. 



Field/Area 


Filter 


N 


expos, time 




Sloan 


T3/MegaP. 


seeing 


comments 


[sq. deg.] 






[s] 


[AB mag] 


Am x 100 


Am x 100 


["] 




WlmOpl 


u* 


5 


3000.51 


25.27 


- 


0.0 ± 0.8 (T) 


1.00 




[wl . -0+1] 












-0.4 ±1.3 (M) 






(0.84) 


g' 


5 


2500.45 


25.55 




-0.4 ± 0.5 (T) 


0.90 
















-7.1 ± 0.9 (M) 








r' 


3 


1500.28 


24.72 




-1.9 ± 1.0 (T) 


0.79 






i' 


7 


4305.67 


24.61 




-0.1 ± 0.8 (T) 


0.85 
















-0.8 ± 0.7 (M) 








z' 


11 


6601.19 


23.88 




0.2 ± 1 .5 (T) 


0.79 
















-7.7 ±1.1 (M) 






Wlm0p2 


u* 


5 


3000.58 


25.35 


- 


-0.6 ± 1.1 (T) 


1.00 




[wl . -0+2] 












O O i 1 C /IV M '\ 

—2.2 ±1.5 (M) 






(0.76) 


g' 


7 


3500.46 


25.72 




-0.5 ±1.0 (T) 


0.95 
















-2.5 ± 0.9 (M) 








r' 


2 


1 1"! 1"! 1"! 1 O 

1 000. 18 


24.61 




-0.5 ± 0.5 (T) 


0.82 


no ch. 21, 35 




i' 


7 


4305.66 


24.72 




1.3 ± 2.4 (T) 


0.74 
















0.2 ±1.8 (M) 








z' 


10 


6000.83 


23.64 




-1.2 ± 3.5 (T) 


0.79 


fr. res. 














-7.3 ±1.5 (M) 






Wlm0p3 


u* 


5 


3000.58 


25.27 


- 


-1.0 ± 1.4 (T) 


0.87 




[wl . -0+3] 












5.3 ±1.2 (M) 






(0.75) 


g' 


5 


2500.49 


25.56 




-0.7 ± 0.8 (T) 


0.90 
















-1.9 ± 0.6 (M) 








r' 


2 


1000.18 


24.62 




-1.2 ± 0.5 (T) 


0.85 


no ch. 21, 35 




i' 


7 


4305.65 


24.59 




-0.5 ± 0.7 (T) 


0.74 
















AO i (\ C /IV T\ 

— U.o ± 0.5 (M) 








„, 
z 


1 f\ 

10 


6000.84 


23.59 




1 O i O T / r T' \ 

-1.3 + 2.7 (1) 


0.82 


fr. res. 














O 1 i 1 A /IV /I^ 

-5.1 ±1.4 (M) 






Wlmlpl 


u* 


7 


4200.60 


25.50 






1.00 




[wl . -1+1] 


g' 


8 


a r\r\r\ i si 

4000.76 


25.73 






0.66 




(0.84) 


r' 


2 


1 f\f\f\ f*\ f\ 

1000.20 


24.51 




- 


0.59 






i' 


9 


5535.73 


24.53 




■ 


0.79 






z' 


6 


3600.40 


23.30 






0.75 




Wlmlp2 


u* 


5 


3000.52 


25.03 




-12.0 ± 0.6 (T) 


0.92 




[wL-1+2] 












-9.2 ± 1.5 (M) 






(0.82) 


g' 


5 


r\ rAA A O 

2500.48 


25.44 




3.3 ± 0.7 (1) 


0.87 
















-4.1 ±1.0 (M) 








r' 


4 


2000.38 


24.80 




1.0 ± 2.2 (T) 


0.87 
















-2.1 ± 2.0 (M) 








i' 


7 


4305.66 


24.52 




-0.1 ± 1.8 (T) 


0.71 
















1 /; i o 1 /iv t\ 

1 .o ± z.l (M) 








z' 


1 1 


6601.18 


23.87 




-1.1 ±1.5 (Ij 


0.71 


rr. res. 














0.3 ± 0.8 (M) 






Wlmlp3 


u* 


5 


3000.50 


24.95 




-18. 3 ± 0.5 (T) 


0.79 




|_W1 . — 1 + 3J 


g 


i 
j 


1 jUU. jU 


2.3. LL 




l.U ± U.o (1 J 


U. / / 




(0.76) 


r' 


2 


1000.23 


24.58 




1.4 ± 0.8 (T) 


0.79 


no ch. 21, 35 




i' 


7 


4305.70 


24.58 




-0.3 ± 1.0 (T) 


0.74 






z' 


10 


6001.00 


23.43 




-11.9 ±2.3 (T) 


0.71 


fr. res. 


Wlplml 


u* 


5 


3000.42 


25.26 


12.2 + 5.5 




0.85 




[wl.+l-l] 


g' 


6 


3000.47 


25.79 


2.7 + 2.8 




0.87 




(0.85) 


r' 


2 


1000.16 


24.50 


4.2 + 3.4 




0.71 






i' 


7 


4305.53 


24.85 


5.2 + 4.8 




0.71 





T. Erben et al.: CARS - Five-band multi-colour data from 37 sq. deg. archival CFHTLS observations 



Field/ Area 


Filter 


N 


expos, time 


m\im 


Sloan 


T3/MegaP. 


seeing 


comments 


[sq. deg.] 






[s] 


[AB mag] 


Am x 100 


Am x 100 


["] 






z' 


6 


3600.40 


23.52 


1.3+5.7 


- 


0.85 




Wlplpl 


u* 


5 


3000.47 


25.28 


- 


0.6 ± 0.9 (T) 


0.92 




[wl.+l+l] 












6.6 ± 0.7 (M) 






(0.75) 


g' 


10 


5000.88 


25.97 


- 


-0.3 ± 0.7 (T) 


0.95 
















-2.0 ± 0.5 (M) 








r' 


2 


1000.18 


24.57 


- 


-1.6 + 0.6 (T) 


0.85 






i' 


7 


4305.72 


24.70 


- 


-0.1 ±0.8 (T) 


0.82 
















-0.5 ± 0.5 (M) 








z' 


12 


7201.21 


23.94 


- 


0.9 + 1.5 (T) 


0.74 


fr. res. 














-9.3 ± 1.3 (M) 






Wlplp2 


u* 


5 


3000.48 


25.29 


- 


-3.9 + 1.0 (T) 


1.00 




[wl. + l+2] 


g' 


5 


2500.42 


25.60 


- 


-1.3 ± 1.0 (T) 


0.87 




(0.83) 


r' 


2 


1000.16 


24.50 


- 


-0.1+0.6 (T) 


0.74 






i' 


6 


3690.57 


24.45 


- 


-0.1 ± 1.0 (T) 


0.77 






z' 


10 


6001.00 


23.60 


- 


-0.4 ± 1.8 (T) 


0.63 


fr. res. 


Wlplp3 


u* 


5 


3000.47 


25.30 


- 


-2.7 + 1.1 (T) 


0.95 


no ch. 21, 35 


[wl.+l+3] 












-1.8 ±1.7 (M) 






(0.80) 


a' 
& 


5 


2500.42 


25.53 




2.4 ± 0.8 (T) 


0.95 


no ch. 21, 35 




r' 


4 


2000.22 


24.72 


- 


-1.4 ± 1.4 (T) 


0.89 


no ch. 21, 35 




i' 


8 


4920.84 


24.63 


- 


-0.4 ± 0.8 (T) 


0.95 


no ch. 21, 35 




z' 


10 


6001.03 


23.63 


- 


-0.9 ± 2.2 (T) 


0.69 


no ch. 21, 35; fr. 














-9.9 ± 1.7 (M) 




res. 


Wlp2pl 


u* 


5 


3000.53 


25.32 


- 


-3.1 ± 1.4 (T) 


1.00 


no ch. 3 1 


[wl.+2 + l] 


g' 


5 


2500.45 


25.56 


- 


-1.9 ± 1.1 (T) 


1.00 


no ch. 31 


(0.84) 


r' 


2 


1000.18 


24.42 


- 


-2.8 ± 0.6 (T) 


0.95 


no ch. 3 1 




i' 


8 


4960.74 


24.63 


- 


0.6 ± 0.6 (T) 


0.90 


no ch. 3 1 




z' 


10 


6000.82 


23.74 


- 


-0.9 ± 2.2 (T) 


0.71 


no ch. 31; fr. res. 


Wlp2p2 


u* 


5 


3000.56 


25.21 


- 


-2.8 ± 1.1 (T) 


0.98 


no ch. 21, 35 


[wl.+2+2] 


g' 


5 


2500.45 


25.63 


- 


1.4 ± 1.2 (T) 


0.95 




(0.81) 


r' 


2 


1000.18 


24.53 


- 


-3.4 ± 0.6 (T) 


0.95 






i' 


8 


4960.76 


24.80 


- 


0.7 ± 1.1 (T) 


0.85 






z' 


10 


6000.95 


23.73 


- 


-1.1 ± 1.3 (T) 


0.74 




Wlp2p3 


u* 


7 


5950.29 


25.61 


- 


-0.9 ± 0.5 (T) 


1.00 


no ch. 3 1 


[wl.+2+3] 


g' 


5 


2500.35 


25.62 


- 


-0.3 ± 0.5 (T) 


0.95 




(0.83) 


r' 


4 


2000.37 


24.82 


- 


-6.3 ± 0.9 (T) 


0.71 






i' 


7 


4340.55 


24.59 


- 


-0.4 ± 0.9 (T) 


0.95 






z' 


9 


7200.41 


23.80 


- 


0.9 ± 2.3 (T) 


0.69 


no ch. 31; fr. res. 


Wlp3mQ 


u* 


5 


3000.55 


25.17 


11.9 + 6.6 


6.4 ± 2.2 (M) 


0.63 


no ch. 3 1 


[W1.+3-8] 


g' 


5 


2500.43 


25.54 


-1.8 ± 3.0 


-2.6 ± 0.6 (T) 


0.87 


no ch. 31; m. ZP 


(0.81) 


r' 


2 


1000.17 


24.33 


2.4 + 3.0 


0.1 ± 1.3 (T) 


0.69 


no ch. 3 1 




i' 


5 


3100.23 


24.44 


3.5 + 3.5 


3.1 ± 1.0 (T) 


0.71 


no ch. 3 1 ; sat. tr. 




z' 


12 


7201.49 


23.59 


-0.4 ± 7.2 


-8.8 ±3.7 (M) 


0.71 


no ch. 31; fr. res. 


Wlp3pl 


u* 


5 


3000.26 


25.31 


- 




0.85 


no ch. 31 


[wl.+3+l] 


g' 


6 


3000.27 


25.62 


- 


0.0 ± 0.5 (T) 


0.95 


no ch. 3 1 


(0.85) 


r' 


2 


1000.10 


24.39 


- 


0.5 ± 0.7 (T) 


0.85 


no ch. 3 1 




i' 


7 


4340.32 


24.60 


- 


-0.8 ±1.1 (T) 


0.95 


no ch. 3 1 




z' 


6 


3600.47 


23.55 


- 




0.69 


no ch. 3 1 


Wlp3p2 


u* 


7 


4200.42 


25.45 


- 


5.0 ± 1.5 (M) 


0.90 


no ch. 3 1 


[wl.+3+2] 


g' 


5 


2500.41 


25.54 


- 


0.3 ± 0.9 (T) 


0.85 


no ch. 3 1 ; m. ZP 


(0.82) 


r' 


3 


1500.22 


24.51 


- 


-1.2 ± 1.2 (T) 


0.83 


no ch. 3 1 




i' 


6 


3720.24 


24.48 




0.6 ± 0.8 (T) 


0.69 


no ch. 3 1 




z' 


6 


3600.58 


23.31 




-7.3 ± 0.9 (M) 


0.55 


no ch. 31 


Wlp3p3 


u* 


4 


2400.39 


25.01 




3.8 ±0.9 (M) 


1.11 


no ch. 3 1 


[wl.+3+3] 


g' 


5 


2500.20 


25.51 




0.1 ±0.8 (T) 


0.95 


no ch. 3 1 



T. Erben et al.: CARS - Five-band multi-colour data from 37 sq. deg. archival CFHTLS observations 



Field/ Area 


Filter 


N 


expos, time 


m\im 


Sloan 


T3/MegaP. 


seeing 


comments 


[sq. deg.] 






[s] 


[AB mag] 


Am x 100 


Am x 100 


["] 




(0.84) 


r' 


2 


1000.10 


24.39 


- 


-2.9 ± 0.6 (T) 


0.87 


no ch. 3 1 




i' 


7 


4340.33 


24.51 


- 


-1.3 + 0.9 (T) 


0.82 


no ch. 3 1 




z' 


6 


3600.62 


23.39 


- 


-7.2 ± 0.9 (M) 


0.55 


no ch. 31 


Wlp4m0 


u* 


5 


3000.54 


25.35 


9.4 + 4.7 


6.4 + 0.7 (M) 


0.72 


no ch. 3 1 


[wl.+4-Q] 


g' 


5 


2500.41 


25.50 


-1.5 + 3.1 


-0.8 ± 1.1 (T) 


0.87 


no ch. 31; m. ZP 


(0.82) 


r' 


2 


1000.19 


24.33 


2.7 + 2.5 


0.3 ± 0.6 (T) 


0.71 


no ch. 3 1 




i' 


4 


2480.15 


24.15 


0.1 + 3.1 


0.9 + 1.1 (T) 


0.95 


no ch. 3 1 




z' 


6 


3600.72 


23.48 


0.4 + 5.9 


-7.1+0.8 (M) 


0.69 


no ch. 3 1 


Wlp4pl 


u* 


5 


3000.54 


25.42 


- 




0.79 




[wl.+4+l] 


g' 


5 


2500.43 


25.45 


- 


-1.5 + 0.7 (T) 


0.85 


no ch. 31; m. ZP 


(0.82) 


r' 


2 


1000.16 


24.21 


- 


1.6 + 0.6 (T) 


0.79 


m. ZP 




i' 


7 


4340.62 


24.48 


- 


0.3 ± 0.8 (T) 


0.79 


m. ZP 




z' 


6 


3600.52 


23.34 


- 




0.63 




Wlp4p2 


u* 


5 


3000.60 


25.29 


- 




0.77 




[wl.+4+2] 


g' 


5 


2500.44 


25.49 


- 


-0.9 ± 1.3 (T) 


0.85 


no ch. 31; m. ZP 


(0.87) 


r' 


2 


1000.16 


24.22 


- 


-2.2 ± 0.8 (T) 


0.74 


m. ZP 




i' 


7 


4340.64 


24.44 


- 


-0.9 + 1.1 (T) 


0.87 


m. ZP 




z' 


6 


3600.50 


23.45 


- 




0.55 




Wlp4p3 


u* 


5 


3000.25 


25.32 


- 


5.9 + 0.8 (M) 


0.98 


no ch. 3 1 


[wl.+4+3] 


g' 


5 


2500.41 


25.65 


- 


-0.9 ± 2.3 (T) 


0.95 


no ch. 3 1 


(0.83) 












-1.4 ± l.l (M) 








r 


2 


1 r\r\r\ 1 h 

1000.17 


A A £ 

24. 46 




A T ■ 1 A /TA 

-0.7 ± 1.0 (1) 


A AO 

0.93 


no ch. 3 1 




• ) 
l 


'"7 

7 


A 1 A n 1 A 

4340.34 


A A A 

24 .44 




1 A i A H / r V\ 

-1.0 ± 0.7 (1) 


A AC 

0.95 


no ch. 3 1 




z 


6 


3600.3o 


23.29 






A -7-7 
0. / / 


no ch. 3 1 


W3mlm2 


u* 


5 


3000.96 


24.78 


30.2 + 3.9 


25.9 + 2.5 (M) 


0.87 




[W3.-1-2] 


g' 


5 


2501.02 


25.48 


0.2 ± 1.9 


0.2 ± 0.6 (T) 


0.87 


m. ZP 


(0.83) 












-4.9 + 1.5 (M) 








r' 


2 


1000.47 


24.51 


-0.1+2.7 


-1.3+0.7 (T) 
-9.9 ± 3.7 (M) 


0.63 






i' 


7 


4306.47 


24.36 


0.6 + 3.8 


-0.6 + 1.3 (T) 


0.66 






z' 


6 


3601.20 


23.41 


-6.1 ±4.2 


-9.4 + 2.3 (M) 


0.67 




W3mlm3 


u* 


5 


3000.98 


24.87 


29.9 + 3.8 


26.9 + 2.6 (M) 


0.74 


no ch. 21 


[W3.-1-3] 


g' 


5 


2500.90 


25.63 


0.4 ± 2.0 


0.1+0.5 (T) 


0.87 


no ch. 21; m. ZP 


(0.82) 












-4.7 + 1.1 (M) 








r' 


2 


1000.28 


24.52 


2.9 ± 2.2 


-2.2+ 1.6 (M) 


0.97 


no ch. 21 




i' 


7 


4306.61 


24.36 


0.6 ± 2.7 


0.3 ± 0.9 (T) 
-4.4+1.6 (M) 


0.66 


no ch. 21; m. ZP 




z' 


6 


3601.17 


23.36 


-3.8+5.6 


-9.1 +2.2 (M) 


0.59 


no ch. 21 


W3m2m2 


u* 


5 


3001.06 


25.33 


29.4 + 4.1 


25.3 + 2.2 (M) 


0.77 




[W3.-2-2] 


g' 


5 


2500.89 


25.59 


0.7 ± 2.2 


-0.3 ± 0.5 (T) 


0.90 




(0.85) 












-3.8 ± 1.0 (M) 








r' 


2 


1000.42 


24.56 


0.7 ± 2.7 


-1.1 ± 1.0 (T) 
-4.8 ± 2.2 (M) 


0.57 






i' 


7 


4306.23 


24.46 


0.3+2.7 


0.2 ± 0.8 (T) 


0.66 






z' 


6 


3601.22 


23.43 


-3.3+4.9 


-8.5 + 1.7 (M) 


0.64 




W3m3m2 


u* 


5 


3000.99 


25.27 


29.0 + 5.0 




0.68 




[W3.-3-2] 


g' 


5 


2500.87 


25.40 


0.4 + 2.2 


-0.2 ± 1.1 (T) 


0.66 




(0.84) 


r' 


2 


1000.40 


24.20 


4.5 ± 2.4 


1.9 + 0.7 (T) 


0.66 






i' 


7 


4306.27 


24.43 


0.1 + 3.3 


0.0 ± 0.9 (T) 


0.53 






z' 


6 


3601.20 


23.30 


-2.3 ± 4.7 




0.55 




W3m3m3 


u* 


5 


3001.02 


25.29 


28.8+4.1 




0.74 




[W3.-3-3] 


g' 


5 


2500.93 


25.52 


0.2 ± 2.3 


-0.5 ± 0.6 (T) 


0.82 




(0.83) 


r' 


2 


1000.42 


24.25 


4.4 ± 2.6 


1.4 + 0.6 (T) 


0.66 






i' 


7 


4306.19 


24.49 


-0.8 ± 3.0 


0.1+0.8 (T) 


0.50 






z' 


6 


3601.21 


23.49 


-3.1+5.5 




0.58 





T. Erben et al.: CARS - Five-band multi-colour data from 37 sq. deg. archival CFHTLS observations 



Field/Area 


Filter 


IN 


expos, time 




Sloan 


1 3/MegaF. 


seeing 


comments 


[sq. deg.] 






[s] 


r a n n 

[AB mag] 


A v * 1 (\(\ 

am x 100 


A v * 1 (\(\ 

am X 100 


[ ] 




W4m0m0 


u* 


5 


3000.26 


25.38 


26.4 + 3.8 


- 


1.03 




[w4.+0+0] 


g' 


5 


2500.40 


25.43 


-0.1 +2.0 


-1.9 ±0.6 (M) 


0.79 




(0.82) 


r' 


3 


1500.22 


24.23 


-0.2 ± 2.2 


1.0 ±0.6 (M) 


0.61 






i' 


7 


4305.65 


24.62 


5.3 + 3.3 


- 


0.71 






z' 


12 


7200.87 


23.78 


-2.3 ± 4.6 


- 


0.66 




W4m0ml 


u* 


5 


3000.29 


25.30 


23.2 + 4.0 


- 


0.74 




[w4.+S-l] 


g' 


10 


5000.68 


25.85 


-0.2 ± 2.5 


-7.7 ± 1.6 (M) 


0.82 




(0.83) 


r' 


2 


1000.10 


24.30 


1.7 + 2.3 


-5.4 ± 1.7 (M) 


0.67 






i' 


7 


4305.36 


24.62 


-0.8 ± 3.0 


- 


0.56 






z' 


6 


3600.27 


23.29 


-3.2 + 5.0 


- 


0.50 




W4m0m2 


u* 


5 


3000.31 


25.41 


24.0 ± 4.5 


- 


0.71 




[w4.+0-2] 


g' 


10 


5000.58 


25.83 


-0.3 ± 1.9 


-6.2 ± 1.2 (M) 


0.77 




(0.81) 


r' 


2 


1000.10 


24.28 


2.1+2.2 


- 


0.58 






i' 


7 


4305.58 


24.77 


0.7 ± 2.6 


- 


0.61 






z' 


6 


3600.34 


23.62 


-4.4 + 4.1 


- 


0.63 




W4mlml 


u* 


5 


3000.24 


25.26 


24.2 + 3.9 


- 


0.69 




[W4.-1-1] 


g' 


13 


6500.96 


26.00 


2.9 ± 2.4 


- 


0.79 




(0.79) 


r' 


2 


1000.16 


24.32 


1.4 + 2.5 


- 


0.87 






i' 


6 


3690.44 


24.62 


-0.2 ± 3.2 


- 


0.71 






z' 


5 


3000.24 


23.28 


-1.2 + 4.8 


- 


0.48 




W4mlm2 


u* 


6 


3600.41 


25.45 


33.2 + 8.1 


- 


0.66 




[w4.-l-2] 


g' 


5 


2500.53 


25.54 


-1.8 ± 1.8 


-7.6 ± 1.3 (M) 


0.79 




(0.82) 


r' 


2 


1000.15 


24.34 


2.9 + 2.4 


- 


0.50 






i' 


7 


4305.43 


24.60 


1.8 + 2.9 


- 


0.72 






z' 


5 


3000.38 


23.37 


-0.3 ± 4.7 


- 


0.51 




W4plm0 


u* 


5 


3000.43 


25.29 


26.2 + 4.1 


- 


0.90 




[w4.+l+8] 


g' 


5 


2500.30 


25.40 


1.6 + 2.1 


-2.3 ± 0.6 (M) 


0.67 




(0.75) 


r' 


2 


1000.16 


24.35 


3.0 + 2.2 


- 


0.94 






i' 


7 


4305.52 


24.31 


0.4 + 3.3 


- 


0.53 






z' 


6 


3600.36 


23.33 


-5.0 + 4.2 


- 


0.55 




W4plml 


u* 


5 


3000.22 


25.33 


27.5 + 3.8 


- 


0.85 




[w4.+l-l] 


g' 


5 


2500.34 


25.33 


1.2 + 2.0 


- 


0.83 




(0.80) 


r' 


2 


1000.17 


24.28 


-0.8 ± 2.5 


- 


0.67 






i' 


14 


8611.03 


24.86 


1.4 + 3.7 


- 


0.66 






z' 


6 


3600.33 


23.35 


-4.2 ± 4.7 


- 


0.63 




W4plm2 


u* 


5 


3000.39 


25.15 


12.7 + 3.9 


- 


0.87 




[w4.+l-2] 


g' 


5 


2500.31 


25.37 


0.9 ± 1.9 


-4.9 ± 1.3 (M) 


0.85 




(0.82) 


r' 


2 


1000.10 


24.34 


1.7 + 2.1 


- 


0.61 






i' 


7 


4305.41 


24.57 


1.1 + 3.4 


-5.8 ± 1.1 (M) 


0.71 






z' 


5 


3000.48 


23.12 


-2.1+5.0 


- 


0.53 




W4p2m0 


u* 


5 


3000.31 


25.24 


28.2 + 4.4 


- 


0.79 




[W4.+2-8] 


g' 


5 


2500.47 


25.30 


-1.6 + 2.4 


-4.4 ± 1.6 (M) 


0.74 




(0.77) 


r' 


2 


1000.14 


24.36 


2.0 + 2.3 


- 


0.63 






i' 


7 


4305.53 


24.68 


0.0 + 3.0 


- 


0.57 






z' 


6 


3600.39 


23.20 


-4.9 ± 3.9 


- 


0.79 




W4p2ml 


u* 


5 


3000.36 


25.18 


19.7 + 3.7 


14.5 ± 2.0 (M) 


0.98 




[w4.+2-l] 


g' 


5 


2500.39 


25.45 


-0.7 ± 1.9 


-6.4 ±1.1 (M) 


0.85 




(0.80) 


r' 


2 


1000.22 


24.12 


1.6 + 2.4 


- 


0.85 






i' 


7 


4305.69 


24.53 


-3.4 + 2.9 


-9.9 ± 0.9 (M) 


0.66 






z' 


12 


7200.72 


23.77 


-8.0 + 4.5 


-15.7 ±2.1 (M) 


0.74 




W4p2m2 


u* 


4 


2400.24 


25.18 


16.4 + 3.7 




1.00 




[W4.+2-2] 


g' 


5 


2500.37 


25.36 


-1.3 + 2.1 




0.77 




(0.83) 


r' 


2 


1000.22 


24.05 


2.5 ± 2.5 




0.90 






i' 


13 


7995.87 


24.98 


1.9 + 4.1 




0.63 





T. Erben et al.: CARS - Five-band multi-colour data from 37 sq. deg. archival CFHTLS observations 



Field/ Area 


Filter N 


expos, time 


m\im 


Sloan 


T3/MegaP. 


seeing 


comments 


[sq. deg.] 




[s] 


[AB mag] 


Am x 100 


Am x 100 


["] 






z' 10 


6000.55 


23.66 


-3.7+4.0 




0.72 





